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SECTION  1 


INTRODUCTION  AND  SUMMARY  OF  RESULTS 


1.1  Statement  of  the  Problem 

Electromagnetic  radiation  can  diffuse  through  a  thin  walled  metallic  enclosure 
and,  perhaps,  interfere  with  the  proper  functioning  of  internal  electronics.  This  report 
shows  how  to  predict  the  current  induced  on  a  wire  within  a  rectangular  metallic  box 
(figure  1)  which  is  irradiated  with  an  el^tromagnetic  pulse.  The  wire  is  connected  between 
opposite  walls  of  the  box  and  the  incident  pulse  has  arbitrary  time  history  and  angle  of 
incidence.  Simple  formulas  are  derived  for  a  straight  wire.  The  case  of  an  arbitrarily 
shaped  wire  within  any  shape  enclosure  is  discussed  semi-quantitatively.  Good  agreement 
has  been  found  between  the  analytic  formulas  and  computer  code  results. 

This  work  was  motivated  by  an  experiment  at  Harry  Diamond  Laboratory  which 
measured  the  current  produced  on  a  wire  within  a  shielded  room  by  an  EMP  simulator. 
The  walls  of  the  room  were  entirely  coated  with  metal  so  that  diffusion  was  the  primary 
mechanism  for  leakage  of  radiation  into  the  enclosure.  A  future  document  will  compare 
the  results  of  these  experiments  with  the  theory  developed  in  this  report.  Comparisons 
are  shown  here  between  the  theory  and  some  code  results. 

It  is  helpful  to  consider  three  distinct  regimes  in  frequency  space^~^.  In  the  low 
frequency  regime  the  skin  depth  is  larger  than  the  wall  thickness  of  the  enclosure.  For 
instance,  if  the  conductivity  and  wall  thickness  of  the  enclosure  are  10^  and  1  millimeter, 
respectively,  then  the  low  frequency  regime  extends  from  zero  frequency  to  about  1.6x10^ 
Hertz  (angular  frequency).  The  intermediate  frequency  regime  extends  up  to  the  frequency 
of  the  lowest  internal  mode  of  the  enclosure.  For  example,  the  lowest  mode  of  a  cube  two 
meters  on  each  side  is  2.5x10*  Hertz.  The  high  frequency  regime  extends  upward  from 
the  lowest  mode. 

Figures  2,  3,  and  4  show  typical  results  yielded  by  the  theory.  Figures  2,  3 
and  4a  show  the  wire  current  as  a  function  of  the  frequency  of  the  incident  pulse  in  the 
low,  intermediate  and  high  frequency  regimes  for  a  box  2  meters  on  each  side  irradiated 
with  an  external  field  of  Hx=l  Ampere/meter.  The  wire  parameters  and  other  details  are 
indicated  in  sections  1.2  -  1.4.  As  expected,  the  wire  current  is  largest  at  low  frequencies 
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because  the  enclosure  is  a  good  shield  if  the  skin  depth  is  smaller  than  the  wall  thickness. 
Figure  4b  shows  the  wire  current  in  the  time  domain.  The  box  is  irradiated  by  a  ten 
nanosecond  long  square  pulse  whose  magnetic  field  rises  abruptly  to  a  peak  Hx=1000 
Amperes/meter.  The  wire  current  rises  to  a  peak  value  of  about  0.5  ma  which  persists 
for  about  100  microseconds.  The  wire  resistance  is  0.01  ohms/meter  and  the  inductance 
is  10~^  henry /meter.  The  other  parameters  are  the  same  as  in  the  calculations  of  figures 
2-4. 


The  remainder  of  this  introduction  (sections  1.2,  1.3  and  1.4)  summarizes  the 
final  results  of  this  analysis.  The  final  analytical  formulas  for  the  three  regimes  are  dis¬ 
cussed  along  with  the  code  comparisons  which  corroborate  them.  Section  1.4  discusses  the 
physics  of  the  wire  response  and  extends  the  analysis  to  arbitrarily  shaped  enclosures  and 
wire  shapes.  Section  1.5  summarizes  some  basic  mathematical  facts  which  are  used  in  sec¬ 
tions  2-5.  Sections  2  through  5  contain  the  detailed  derivations  of  the  results  sununarized 
in  section  1.  Section  6  discusses  the  code  HFC  which  treats  external  scattering  from  a  box. 
A  previous  report  discusses  the  low  frequency  regime  in  greater  detail.  For  the  reader’s 
convenience  it  is  reproduced  here  (with  some  minor  corrections)  in  Appendix  A. 

MKS  units  are  used  in  this  report. 


1.2  Low  Frequency  Regime 

This  subsection  summarizes  the  final  results  of  the  low  frequency  regime.  The 
detailed  derivation  are  shown  in  Appendix  A. 

For  boxes  of  any  practical  size  the  wavelength  of  the  incident  electromagnetic 
wave  in  the  low  frequency  regime  is  much  larger  them  the  box.  Consequently,  in  the 
neighborhood  of  the  box  the  incident  magnetic  field  is  spatially  uniform  and  is  given  by 
the  real  part  of 


H°e-‘"‘  .  (1-1) 

Here,  a  single  Fourier  component  of  the  incident  pulse  will  be  considered  which  oscillates 
at  the  frequency  u;.  The  time  domain  wire  current  will  be  considered  below. 
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The  uniform  external  field  can  be  resolved  into  its  vector  components  H° , 
and  along  the  three  axes  in  figure  1.  Then  the  current  in  the  wire  is  given  by  the  final 
result 


IM  = 


,  C,  b 
7-iwL^°  2^b  +  c 


iu;MoH;  a  b 
7-iwL^^  2^a  +  b 


(1-2) 


Here,  the  coordinates  of  the  wire  are  x  =  xq  and  z  =  Zq.  H*  and  H*  are  the  vector 
components  of  the  magnetic  field  inside  the  box 


h; 


1  —  iw/wi 


1  —  lUi/hJi 


(1-3) 


which  are  also  spatially  uniform  (Appendix  A  and  section  2).  Also  a,b,c  are  the  box 
dimensions  and 


The  resistance  per  unit  length  of  the  wire  (at  zero  frequency)  is,  of  course®. 


Ro  = 


1  .  ohm  ^ 

JTWoao  'meter' 


(1-4) 


in  terms  of  its  radius  Wq  and  conductivity  ctq.  But  at  higher  frequency  where  the  skin 
depth 


tfo  = 


i/ 


2 

HoUiOo 


is  less  than  the  wire  radius  (fio  <  wq)  the  wire  impedance®  is 


(1-5) 
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Figure  1.  Problem  geometry. 
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(:-6) 


„  /.  Wo ,  ohm 


The  symbol  7  is  defined  to  be 


_  /  Ro  if  <  Wo 

i  Ro(l-0s^  if^o>wo 

L  is  the  inductance  per  unit  length  of  the  wire  in  the  box.  Above  the  frequency 


(1-7) 


= 


HoffAS. 

where  a  and  A  are  the  wall  conductivity  and  thickness  the  inductance  is 


(1-8) 


L  = 


8mo 

acx* 


00  Sin  ( - )s>n  ( - ) 

^  ^  g  ' 

(;)*  +  (?)* 


.henry. 
'  meter  ^ 


(1-9) 


Below  this  frequency  it  is  approximately  the  inductance  which  the  wire  would  have  in  free 
space  outside  the  box.  The  symbols  Mo  ^nd  Cq  have  their  usual  meaning. 

Figure  2  compares  the  wire  current  predicted  by  equation  1-2  with  the  current 
computed  by  the  computer  code  B0X4  discussed  in  Appendix  A.  The  agreement  is  within 
a  factor  of  four  over  five  orders  of  magnitude  of  frequency.  In  this  example  the  box  is  two 
meters  on  each  side  and  the  walls  are  one  millimeter  thick.  The  wall  conductivity  is  10^ 
mho/meter.  The  wire  resistance  Ro  =  0.0183  ohm/meter  and  it  is  located  one-half  meter 
from  the  x  =  0  and  z  =  0  walls.  The  wire  is  one  centimeter  in  radius  and  the  external 
magnetic  field  is  pointing  in  the  x  -  direction  and  has  a  magnitude  on  one  ampere/meter. 

Table  I  indicates  that  the  wire  current  is  predicted  at  different  positions  of  the 
wire  at  one  hertz  angular  frequency.  Agreement  is  within  twenty  percent  when  the  wire 
location  is  varied  from  Zq  =  0.05  to  Zq  =  0.95  meters.  All  other  parameters  are  identical 
to  the  figure  2  case. 
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Table  1.  Wire  current  as  a  function  of  wire  position. 
Wire  Position  Wire  Current  Wire  Current 


X  (cm) _ Theory  (Amp)  Code  (Amp) 


5 

3.3  xlO-* 

3.5  XIO-® 

10 

3.1  X10-® 

3.3  xlO-® 

50 

1.7  xlO-' 

1.6  xlO-® 

95 

1.7  xlO-® 

1.4  xl0-« 

Table  2  indicates  good  agreement  between  the  analytic  expressiorjs  and  the  code 
when  the  frequency  is  10^  hertz.  The  analysis  always  underpredicts  the  current  by  about 
of  factor  of  three. 


Table  2.  Wire  current  as 


Wire  Position 

Wire  Current 

X  (cm) 

Theory  (Amp) 

5 

0.033 

10 

0.027 

50 

0.011 

95 

0.0011 

a  function  of  wire  position. 


Wire  Current 

Ratio 

Code  (Amp) 

Code  -r  Theory 

0.10 

3.0 

0.081 

3.0 

0.043 

3.9 

0.0029 

2.6 

Figure  2  illustrates  the  general  behavio'  of  the  wire  current.  It  increases  linearly 
with  frequency  to  a  peak  value.  At  sufficiently  high  frequency  the  inductive  reactance  of 
the  wire  becomes  large  enough  to  cause  the  current  to  decrease. 

The  wire  current  can  be  obtained  as  a  function  of  time  by  taking  the  inverse 
fourier  transform  of  equation  1>2.  If  the  incident  pulse  has  an  arbitrary  time  history  given 
by  the  function  then  the  wire  current  is 

I(t)  =  — rdt'H^l^e- 

Ro  —  WjL  70  I L  J 
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I  Wife  Current  (Amperes) 


Ro  “■ jo  IL  J 

In  general,  this  integral  smoothes  out  the  fluctuations  in  the  incident  pulse  so 
that  the  wire  current  is  much  l^s  rapidly  varying.  When  the  incident  pulse  first  reaches 
the  box  the  wire  current  slowly  builds  and  after  the  pulse  has  past  the  box  the  wire  current 
slowly  decays  to  zero. 

These  results  are  explained  in  section  2  and  Appendix  A. 


1.3  Intermediate  Frequency 

The  analytical  formulas  for  the  wire  current  in  the  intermediate  frequency  regime 
are  summarized  here.  The  details  are  in  section  2. 

In  the  intermediate  frequency  regime  equation  1-3  applies  but  the  internal  mag¬ 
netic  fields  within  the  box  are  no  longer  given  by  equation  1-3.  Instead,  the  appropriate 
expressions  to  use  arc 


K 


Eg 

cos(^A)  —  i/?aBin(/?A) 


h: 


HP 

cos(/?A)  —  |/?5sin(/?A) 


(1-11) 


The  wail  thickness  is  A  and  the  symbol  f3  is  the  attenuation  constant  of  a  wave  propagating 
in  the  metallic  wall.  In  terms  of  the  skin  depth  S  in  the  wall  it  is  the  complex  number 


0  = 


1  +  i 
6 


(1  -  12) 


Figure  3  shows  the  intermediate  frequency  regime  wire  current  for  the  case  illus¬ 
trated  in  figure  2.  Clearly,  at  high  frequency  the  wire  current  becomes  very  small  because 
the  walls  are  many  skin  depths  thick  and  the  enclosure  is  a  good  shield. 
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The  pulse  starts  to  diffuse  into  the  enclosure  in  a  time  t^  =  seconds  which 

is  frequently  many  microseconds.  If  the  incident  pulsewidth  is  short  compared  to  this  time 
then  the  wire  current  in  the  time  domain  is  approximately 


i(t)  = 


[i  _  IM  4. 

\tA/  L  WaA  \/5r£itA  \  t  / 


i/a 


{(»o  -  §)r^a>,  -  (xo  -  *)-^H 


b  4-  c 


2'a  +  b 


-1 


t  <  tA/10 


where  «.  = 


Ro  ■"  uJiLi  I  L  J 


*»Wio 


The  wire  current  does  not  depend  upon  the  details  of  the  time  history  of  the  incident  pulse. 
Instead  only  the  time  integrated  imptilse 


Hex  =  r  Hjdt 
00 

Htt.  =  rH?dt 

y-oo 

is  important.  The  above  early  time  approximation  is  valid  provided  Rot^sh/IOL  <  1.  In 
general,  the  wire  current  slowly  rises  to  a  peak  at  about  the  time  r^/lO.  It  then  falls  to 
zero.  These  expressions  make  iise  of  the  early  and  late  time  approximations  for  the  internal 
magnetic  field  which  was  derived  by  Bedrosian  and  Lee’*’. 
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1.4  High  Frequency  Regime 


The  analytic  formulas  for  the  high  frequency  wire  current  are  summarized  here. 

If  a  plane  wave  is  normally  incident  on  the  z  =  c  wall  of  the  box  (figure  1)  then 
the  wire  current  is 


1  = 


-8Z»HL  ^  fni 

^  t,n»l  I  t 


(cosn?r)(l  -  cos 


—  cos  tx)  1  .  tr 
— - r  >  sin  — 

J  ■ 

sin*  gjjji  am 


/ttxo  .  nTTZo 

sin - sm - 

a  c 


7- 


4ia> 


Here  the  Incident  wave  magnetic  field  (at  z  =  c) 


(1  ~  13) 


(1  -  14) 


is  in  the  x  «  direction  and  the  transfer  impedance  of  the  wail  is 


Zt  = 


/? 

<rsin()9A) 


(1  -  15) 


The  speed  of  light  in  vacuum  is  denoted  by  v  and  k  =  u;/v  is  the  wavenumber. 

Figure  4  graphs  the  high  frequency  wire  current  given  by  equation  1-13  for  the 
box  referred  to  in  the  figures  2  and  3  example.  As  expected  the  wire  current  is  minute. 
In  the  high  frequency  regime  the  incident  frequency  can  equal  the  frequency  of  one  of  the 
internal  modes  of  the  box.  One  might  expect  that  this  would  lead  to  a  resonance  which 
might  cause  a  dramatic  change  in  the  wire  current  at  the  resonance  frequency.  Figure  4 
shows  one  such  resonance  near  2  gigahertz  angular  frequency. 

At  low  and  intermediate  frequencies  the  above  formulas  apply  for  any  angle  of 
incidence.  However,  the  case  of  arbitrary  angle  of  incidence  at  high  frequency  is  much  more 
complicated  to  treat  and  no  simple  formula  applies.  Instead  a  computational  method  has 
been  derived.  The  code  HFG  (section  6)  calculates  the  current  on  the  outside  of  the  box. 
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Wire  Current  1  ( 
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These  results  are  then  used  by  the  method  described  in  sections  4  and  5  to  compute  the 
wire  current. 

As  a  practical  matter  the  wire  current  is  so  small  at  high  frequencies  that  a  very 
prediction  is  probably  not  necessary.  An  approximate  method  which  is  convenient 
inploy  is  for  most  applications  more  useful.  It  is  recommended  that  equation  1-13  be 
used  to  yield  an  approximation  in  the  case  that  the  frequency  is  so  high  that  the  wavelength 
of  the  incident  wave  is  much  less  than  the  box  dimensions.  At  lower  frequencies  (in  the 
high  frequency  regime)  the  wire  current  might  be  estimated  with 


-i^vH*jSin(ki/) 


4eiM 


cos  ka  sin  ka 


(1  -  18) 


Equation  1-16  furnishes  the  wire  current  for  an  enclosure  formed  by  two  parallel 
plates  of  infinite  extent  (figure  5).  A  pulse  given  by  equation  1-14  is  incident  normally  on 
each  plane.  The  planes  are  a  distance  2a  apart  and  the  wire  coordinate  is  i/q  (figure  5). 

Figure  5  also  graphs  the  wire  current  for  a  parallel  plate  enclosure  formed  by 
plates  2  meters  apart  having  the  same  conductivity  and  thickness  as  in  the  previous  exam¬ 
ple.  A  similar  wire  is  located  0.5  meters  from  one  plane.  The  total  field  on  the  planes  is 
one  ampere/meter.  The  magnitude  of  the  current  is  about  the  same  as  in  the  box  example 
in  figure  4  but  the  resonance  spikes  are  more  visible. 

The  code  HFC  has  been  written  to  treat  the  high  frequency  case  with  a  great 
deal  of  precision.  Figure  6  shows  that  it  calculates  correctly  the  scattering  of  an  incident 
electromagnetic  wave.  The  surface  current  density  (amperes/meter)  on  the  outside  of 
the  box  is  graphed  around  the  perimeter  of  the  box.  Agreement  is  seen  to  be  good  with 
published  results  for  the  two  frequencies  discussed  in  section  6.  The  k^  =  1  case  corresponds 
to  a  wavelength  equal  to  three  times  the  box  dimension  a. 

The  wire  current  can  be  computed  as  a  function  of  time  for  the  intermediate  and 
high  frequency  cases  by  taking  the  inverse  fourier  transform 


I(t)  =  ^ 
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of  equations  1>2,  1-13  or  1-16. 

The  parallel  plate  enclosure  is  discussed  in  section  3.  Other  high  frequency 
results  are  explained  in  sections  4,  5  and  6. 


1.5  Physics  of  Wire  Current 

The  fundamental  physics  of  the  wire  coupling  is  easily  understood.  The  metallic 
conductor  very  effectively  shields  the  incident  electric  field  from  the  interior  of  the  con¬ 
ductor  at  all  frequencies.  At  low  frequencies  charge  appears  on  the  external  surface  (the 
conductor  polarizes)  which  creates  an  electric  field  which  cancels  the  incident  field  within 
the  enclosure. 

The  incident  magnetic  field,  however,  slowly  diffuses  through  the  walls.  The  curl 
of  the  magnetic  field  creates  a  small  electric  field  which  drives  current  in  the  metal.  The 
current  creates  its  own  magnetic  field  which  cancels  the  incident  field  so  that  the  internal 
magnetic  field  is  small.  This  is  the  mechanism  of  inductive  shielding  which  occurs  in  the 
low  frequency  regime.  At  higher  frequencies  the  internal  magnetic  field  is  further  reduced 
as  the  wave  propagates  many  skin  depths  through  the  metal. 

The  dotted  line  in  figure  7  defines  a  closed  loop  which  runs  through  the  wire 
and  the  metallic  walls.  Around  this  loop  Faraday’s  law  states 

EMF  =  ^  E  •  df  =  iw/io  j  H  •  ndS 

This  equation  indicates  that  the  electromotive  force  driving  current  through  the  wire  equals 
the  rate  of  change  of  magnetic  fiux  through  the  loop.  Clearly,  the  magnetic  flux  increases  as 
the  magnetic  field  increases.  The  equation  also  predicts  that  as  the  frequency  u  increases 
the  EMF  and  wire  current  will  increase.  The  rising  portion  of  the  wire  current  in  figure  2 
is  due  to  this  effect. 

However,  this  source  of  EMF  drives  current  against  two  types  of  impedance. 
Firstly,  the  frictional  force  exerted  by  the  metallic  lattice  impedes  the  electrons  fiow- 
ing  through  it.  This  is  quantified  by  the  wire  resistance.  Secondly,  the  current  fiowing 
through  the  wire  circuit  creates  a  magnetic  field  which  reduces  the  flux  in  the  loop  and, 
consequently,  the  EMF.  This  can  be  quantified  in  terms  of  the  wire’s  self-inductance  L.  As 
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the  freqtiency  increases  the  inductive  reactance  increases  and  so  the  wire  current  decrezises. 
This  causes  the  drop  in  current  in  figure  2  at  the  higher  frequencies. 

In  the  low  frequency  regime  the  internal  magnetic  field  is  largely  independent  of 
the  shape  of  the  enclosure  for  reasons  discussed  in  Appendix  A  and  section  2.  Furthermore, 
the  inductance  L  is  always  of  the  order  of  henry/meter  regardless  of  whether  the  wire 
is  straight  or  (within  reason)  curved.  Consequently,  the  EMF  in  the  loop  in  figure  7  and 
the  wire  current  itself  is  expected  to  be  somewhat  insensitive  to  the  precise  shape  of  the 
enclosure  or  the  wire. 

The  nature  of  the  internal  magnetic  field  in  the  low  frequency  inductive  regime 
has  been  very  well  developed  elsewhere^’^'*^ .  If  a  short  pulse  is  incident  on  an  enclosure  of 
any  shape  the  internal  magnetic  field  slowly  rises  from  zero  and  reaches  a  peak  after  about 
Ti^jA  where 


seconds 

Here  A  is  the  wall  thickness.  This  is  the  time  for  the  radiation  to  diffuse  through  the  walls 
and  is  often  tens  of  microseconds.  The  field  then  decays  to  zero  exponentially  with  a  fall 
time  of 


Td  ~ 


where  f  is  the  ratio  of  the  volume  of  the  enclosure  to  the  volume  of  the  walls.  In  general, 
then,  if  the  incident  pulse  is  shorter  than  then  the  internal  magnetic  field  slowly  rises 
to  a  peak  after  r^/A  seconds.  It  then  falls  even  more  slowly.  The  internal  field  behavior  is 
totally  insensitive  to  the  detailed  time  history  of  the  incident  pulse.  The  internal  electric 
field  is  induced  b  the  changing  internal  magnetic  field  and,  therefore,  is  proportional  to 
the  rate  of  change  of  the  magnetic  field. 

The  wire  current  can  be  expected  to  rise  to  a  peak  at  about  the  time  rp  and 
then  fall  to  zero.  This  behavior  obtains  as  long  as  the  pulsewidth  is  shorter  than  and 
is  independent  of  the  details  of  the  incident  pulse. 


1.6  Mathematical  Remarks 

Several  mathematical  facts  which  will  be  used  below  are  discussed  in  this  section. 
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The  electric  and  magnetic  fields^’®  can  be  written  in  terms  of  the  vector  A  and 
scalar  (f>  potentials 


B  =  V  X  A 

_  ,  dA 

E  =  ■ 

In  this  report  the  Lorentz  gauge  is  used  exclusively. 

Consider  a  scalar  function  0  satisfying  the  differential  equation 


1 

v2  dt^ 


=  0  . 


(1-17) 


(1  -  18) 


Then,  the  transformation 


Ao  =  A  —  Vijt 
drb 

*  = 

yields  new  potentials  which  lead  to  the  same  E  and  B  (equation  1-1)  and  satisfy  the  Lorentz 
gauge  conditions.  Although  A  and  0  must  satisfy  boundary  conditions  on  the  box,  0  can 
have  any  value  on  the  boundaries  (as  well  as  any  initial  conditions)  because  the  electric 
and  magnetic  fields  on  the  boundary  given  by  equation  1-1  are  completely  independent 
of  0.  Therefore,  0  can  be  any  solution  to  the  homogeneous  wave  equation  (equation  1-2) 
which  will  be  seen  to  be  very  useful. 

The  freedom  afforded  by  0  allows  one  to  choose  a  scalar  potential  0  which  is  zero 
on  the  walls.  For  if  0  were  not  zero  on  the  walls  then  0  could  be  chosen  to  have  the  value 

0  =7  —  on  the  walls.  Then  0  must  be  that  solution  of  the  wave  equation  (equation  1-18) 
iw 

which  has  the  prescribed  value  on  the  wall  which  is 
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0  =  /  $.VG  •  dS  ,  (1  -  19) 

J  \U) 

Here  G  is  the  Green’s  function  of  the  wave  equation  satisfying  Dirichlet  boundary  condi¬ 
tions  and  the  integral  is  taken  over  the  wall  surfaces.  Clearly,  another  scalar  potential  i^o 
could  be  defined 


^0  =  ^  + 


dy() 

dt 


Of  course,  i^o  =  0  on  the  walls. 

The  Helmholtz  equation 


(J.  -  20) 


V  V  +  =  —Anp 

arises  from  fourier  transform  of  the  wave  equation  as  well  as  the  diffusion  equation.  It  has 
the  general  solution^ 

0(f)  =  I  pGdV  +  ^  / [GV0  -  0VG1  •  dS 
where  the  Green’s  function  satisfies 

V^G +k^G  = -47r6(f-f')  . 


It  is 


where  0n  S'i'o  eigenfunctions  of  the  Helmholtz  operator 


VV.+kV.  =  0  . 


They  are  normalized  according  to 

/  0m^nd?  =  <5iiin 

0,  G  and  0„  all  satisfy  the  same  boundary  conditions. 
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SECTION  2 


INTERMEDIATE  FREQUENCY  REGIME 


2.1  Introduction 


In  the  intermediate  frequency  regime  the  frequency  is  sufficiently  high  that  the 
skin  depth  is  smaller  than  the  wall  thickness  but  sufficiently  low  that  the  lowest  mode  of 
the  box  is  not  appreciably  excited.  For  instance,  if  the  wall  conductivity  and  thickness  are 
10^  mho/meter  and  one  millimeter,  respectively,  and  the  box  is  2  meters  on  a  side  then 
the  intermediate  frequency  regime  extends  from  about  1.6  xlO^  to  2.5  xlO^  sec'^angular 
frequency. 

At  any  point  within  the  wire  the  current  density  J  is  related  to  the  total  electric 
field  Et  by 


J  =  oqEt  (2  —  l) 

where  Cq  is  the  wire  conductivity.  At  low  frequencies  the  current  is  uniform  throughout 
the  wire  of  radius  Wo  so  the  total  current  is 


I  =  ttwqJ 


(2-2) 


and 


IRo  =  Et 


where  the  wire  resistance  per  unit  lengtlT  is 


Ro- 


1 

TTwgoTo 


ohms 

m 


(2-3) 


At  higher  frequencies,  however,  the  current  runs  primarily  within  a  skin  depth  ^  of  the 
wire  surface  and 
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(2-4) 


!Ro(l  -  0^  =  Et 


Th«  total  electric  field  is 


Et-E*  +  E’'  .  (2-5) 

Here,  E*  is  the  field  within  the  box  when  the  wire  is  not  present  but  the  box  is  irradiated 
with  the  external  electromagnetic  field.  Similarly,  E*'  is  the  field  within  the  box  caused  by 
a  current  I  which  runs  down  the  wire  and  around  the  box  when  the  external  driving  field 
is  absent.  Equations  2-4  and  2-5  indicate  that  the  current  on  the  wire  must  satisfy 


IRo(l-i)^  =  E*  +  E'"  (2-6) 

where  the  fields  are  evaluated  at  the  wire.  The  derivation  of  this  equation  is  very  similar 
to  the  approach  used  in  many  scattering  problems  and  is  explained  in  greater  detail  in 
section  A2  of  Appendix  A. 

E*  and  E*'  are  calculated  in  the  following  two  subsections.  Then  equation  2-6 
will  be  employed  to  calculate  the  current  I  in  the  wire. 


2.2  Calculation  of  E*' 


Figure  7  illustrates  the  problem  geometry.  A  wire  carrying  a  current  density  J 
connects  opposite  walls  of  a  metallic  box.  The  electric  field  within  the  box  is  sought.  Of 
course,  in  order  to  drive  the  current  there  mrat  be  a  source  of  electric  field.  For  present 
purposes,  however,  it  is  only  necessary  to  calculate  the  fields  produced  by  the  current  so 
the  details  of  this  source  are  irrelevant. 

■V 

The  current  in  the  wire  can  be  written  as  a  cosine  series 


(2-7) 
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where  Jp  is  the  Fourier  coeiHcient  and  is  a  function  of  the  radial  position.  The  time 
dependance  is  exp  (-i  uj  t).  The  cosine  series  will  faitnfully  converge  to  any  possible  wire 
current  in  the  open  interval  (0,b)  but  at  the  wire  ends  the  series  fails.  This  causes  no 
difRculty. 

The  Hertz  Vector 

It  is  helpful  to  employ  the  Hertz  vector  FI.  The  polarization  vector  P  is  defined 
in  terms  of  the  current  dj'nsity  by 


di 


-  J 


(2-8) 


By  the  continuity  equation,  therefore,  the  charge  density  on  the  wire  must  be 


/>=  -v  .p  . 


Maxwell’s  equations  then  become 


V  X  H  = 


d%  ^  di 


V.E=  --VP 
Co 


V.B  =  0 


V  X  E  =  - 


£8 

di 


(2-9) 


(2  -  10) 


It  is  easily  shown  that 


these  are  satisfied  by  the  electric  and  magnetic  fields  given  by 


e  =  v(v  •  n)  - 


1  d^ti 

V*  dt* 


B  = 


(2  -  11) 


where  the  Hertz  vector  satisfies 
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(2  -  12) 


V»ft  - 

V*  dt* 


The  speed  of  light  in  vacuum  is  called  v  since  c  denotes  the  box  length. 
In  terms  of  the  vector  and  scalar  potential  (section  A) 


E  = 


(2  -  13) 


B=  7  X  A 


(2  -  14) 


Comparing  equations  2-11  and  2-13  and  2-14  shows  how  11  depends  upon  A  and  <f> 


^  1  dn 

V*  dt 

<t>  =  -V-n  (2-15) 


Boundary  Conditions 

Since  the  wail  conductivity  is  finite  the  inside  and  outside  of  the  box  are  coupied. 
An  exact  solution  would  solve  for  the  exterior  and  interior  fields  together.  Because  the 
skin  depth  in  the  intermediate  region  is  small,  however,  the  coupling  is  extremeiy  small 
and  the  exterior  fields  are  minuscule  compared  to  those  in  the  interior.  Consequently,  the 
walls  will  be  approximated  as  perfect  conductors  to  calculate 

The  Hertz  vector  has  three  components  but  only  the  Ily  component  is  non-zero. 
This  appears  reasonable  since  P  has  only  a  y-component  so  that  only  Ily  is  driven  through 
equation  2-12  But  this  is  insuifirient  to  prove  that  IIx  and  H,  are  non-zero  because  the 
homogeneous  equation  has  non-zero  solutions.  Instead,  the  assumption  that  fix  and  IT, 
are  zero  can  be  justified  at  the  end  of  the  calculation  by  verifying  that  the  final  solution 
satisfies  Maxwell’s  equations.  It  has  been  pointed  out  in  section  A  that  it  is  always  possible 
to  choose  a  scalar  potential  ^  such  that  on  the  walls 
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^  =  0  .  (2  -  16) 

On  the  x=0  and  x=a  and  z=0  and  z=c  walls,  Ey  —  0.  From  equation  2-11  and 
equation  2<15  on  the  walla 


vn  =  0 


(2  -  IT) 


and 


E,  =  ^(r.n)  + 


(2  -  18) 


SO  that  on  these  walls 


n  =  0  .  (2-19) 

On  the  ys=0  and  ysb  walls  equation  2-15  shows 

V  •  n  =  0 


so  that 


an 

<9y 


=  0 


(2  -  20) 


Equations  2-19  and  2-20  are  the  boundary  conditions  on  11. 

Therefore,  by  the  remarks  in  the  last  subsection  of  section  1  the  Hertz  potential 
is 


(2  -  21) 


26 


where  G  is  the  Green’s  function 


00 


G  =  47r 


abc 


sin  ^  cos 


sin  sin 

s. 


moO 


(kS  -  k») 


coe  sjpi  sin 


-a)'- 


(2  -  22) 


where 


kj  =  (7)’  +  (“)’  +  (f )’ 

and  k  =  (j/c.  The  exponent  ^  is  one  if  m=0  and  zero  otherwise.  Here  G  satisfies  Dirichlet 
boundary  conditions  on  two  walls  and  Neumann  conditions  on  the  other.  Of  course  the 
current  runs  primarily  within  a  skin  depth  of  the  surface  of  the  wire  but  the  analysis  is 
simplified  by  approximating  the  current  to  run  along  a  very  thin  filament  so  that 


Jy  =  H  •  (2-23) 

pa»l  “ 

From  equation  2*8  the  polarization  Py  is  found  and  the  Hertz  vector  is 


-4  ^  sin(^)  coa(ggg)  sin(np)  sin(^)  sin(afa) 
ihi  ^  k»  -  k> 

masO 


(2  -  24) 


From  equation  2-11  the  electric  field  within  the  cavity  can  be  computed.  At  the  wire 
location  it  is 


E*(y)  =  t  I ij'^  *7/' .in’  ^  sin'  HI  | 

iwacco  k«  -  k’  a  cl 

msO 


COS 


mzy 


(2  -  25) 


This  is  the  desired  field  at  the  wire  which  will  be  needed  below. 
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The  advantage  of  the  Hertz  vector  approach  is  that  it  simplifies  the  boundary 
conditions.  The  boundary  condition  which  was  employed  is  that  the  electric  field  tangent 
to  the  walls  is  zero  while  the  derivative  of  the  normal  electric  field  is  zero  (since  the 
divergence  of  the  electric  field  is  zero  just  inside  the  wall  where  the  tangential  electric  field 
vanishes).  But  the  electric  field  is  due  to  the  combined  affect  of  both  A  and  Hence, 
the  boundary  condition  mixes  the  behavior  of  A  and  <f>.  The  Hertz  vector  provides  one 
approach  to  unmix  the  boundary  conditions.  Another  possible  approach  is  to  solve  the 
scalar  retarded  potential  equation  subject  to  the  conditions  that  <t>  is  zero  on  the  walls. 
With  this  value  of  the  boundary  condition  for  A  can  be  found  from 


aA 

dt 


-E-W  • 


2.3  External  Driver  E* 


In  this  section  the  electric  field  E*  produced  within  the  enclosure  is  calculated 
when  the  wire  is  not  present.  In  sections  4  and  5  this  problem  is  solved  with  the  help 
of  the  code  HFC  very  precisely  for  the  high  frequency  regime.  This  numerical  method  is 
applicable  to  the  intermediate  frequency  case  also  and  will  yield  a  very  accurate  solution. 
However,  in  the  intermediate  frequency  regime  the  code  treatment  is  unnecessarily  com¬ 
plicated.  The  method  developed  in  this  section  is  expected  to  yield  a  simple,  accurate 
result. 


Figure  8  illustrates  three  enclosures  which  are  irradiated  by  an  electromagnetic 
wave  whose  wavelength  in  free  space  is  much  smaller  than  the  enclosure.  Since  the  wave¬ 
length  is  large  the  external  magnetic  field  is  constant  around  the  enclosure  (Appendix  A) 
and  is  given  by 


RDe-iwt 


(2  -  26) 


The  internal  field  H'  is  constant  in  each  case  and  points  in  the  direction  of  the 
external  field.  For  the  sphere,  parallel  plates,  and  two  polarizations  of  the  cylinder^*’’ 


H‘  =  H° 


(2-27) 
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Figure  8.  (a)  Two  parallel  platee,  (b)  a  cylindrical  fhell  with  longitudinal 

(c)  a  cylindrical  chell  with  transverse  H°,  and  (d)  a  spherical  shell. 
All  cavity  walls  have  thickness  conductivity  a  (taken  from  EMF 
Interaction:  K.  S.  Lee,  AFWIfTR-80-403). 
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H‘  ==  {coa{0A)  -  0Kdn{0A)}'^ 


(2-28) 


H‘  =  H‘* 


(2-29) 


where  0^  =  luj/ioa  {0  =  (1  +  i)/6  where  S  is  the  skin  depth).  The  dimension  S  is  defined 
in  the  figure  and  A  is  the  wall  thickness.  Clearly,  the  internal  field  is  almost  the  same 
regardless  of  the  geometry.  It  is  reasonable  to  assume,  therefore,  that  the  field  within  a  box 
is  also  approximately  uniform  and  its  magnitude  is  roughly  that  of  a  spherical  enclosure. 
If  the  sides  of  the  box  are  equal  the  box  might  best  be  replaced  by  a  sphere  whereas  if  one 
side  of  the  box  is  very  small  it  would  best  resemble  two  parallel  plates.  Fortunately,  the 
above  formulas  indicate  that  all  geometries  have  about  the  same  internal  fields. 

Equation  2-27,  therefore,  will  be  used  to  approximate  the  uniform  field  within 
the  box.  The  sphere’s  radius  is  chosen  to  be 


(2-30) 

SO  that  the  volume  of  the  sphere  and  box  are  the  same.  In  Appendix  A  the  field  within  a 
box  and  sphere  are  shown  to  be  nearly  the  same  in  the  low  frequency  regime  by  numerically 
solving  for  the  internal  field  with  the  code  BOX4. 

Equation  2-27  is  derived  in  the  remainder  of  this  subsection.  In  section  2  of 
Appendix  A  it  is  shown  that  the  magnetic  field  within  an  enclosure  is  governed  by  the 
static  magnetic  field  equations  whenever  the  enclosure  is  driven  at  a  frequency  much  less 
than  the  lowest  modal  frequency  of  the  box.  Consequently  the  internal  (and  external) 
magnetic  field  are  given  by  the  gradient  of  a  magnetic  potentiai  which  satisfies  Laplace’s 
equation.  Assuming  that  the  form  of  the  potentials  are  (c.f.  figure  9) 


n‘"  =  ffreostf 

f^oui  ^  H°rcos+^cos^  (2-31) 

the  fields  tangential  to  the  sphere  are 
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H4  =  -ffsintf 

Hi  =  sin^  —  ^  sind  .  (2-32) 

2t 

Within  the  conductor 


V*H  +  /?*H  =  0  (2  -  33) 

which  can  be  solved  by  separation  of  variables  to  be  the  product  of  spherical  Hankel 
functions  and  sin  6.  In  the  limit  that  the  skin  depth  is  smaller  than  the  spheres  radius  the 
two  independent  solutions  are 

sin(d)  and  —  sin(tf)  (2  —  34) 

The  electric  field  associated  with  each  magnetic  field  is  given  by 

E  =  -V  X  H  (2-35) 

a 

so  that  the  two  independent  solutions  of  electric  field  are 

8in(g)  and  ^5ie“‘^'sin(^)  (2  —  36) 

ax  ax 

The  continuity  of  the  magnetic  and  electric  field  across  the  sphere’s  inner  surface 

yields 


i3s. 


(2  -  37) 


aS.  ai 


(2  -  38) 
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whereas  the  continuity  of  the  magnetic  field  across  the  outer  surface  provides 


5*  y3b  I3h 


(2  -  39) 


Within  the  conductor  the  normal  component  of  the  magnetic  field  satisfies  Fara¬ 
day’s  law.  At  the  outer  surface  it  equals  the  normal  component  of  the  external  magnetic 
field.  Consequently, 


aoout 

(V  X  E)„  =  iw/ioHn  =  luifio  (2  -  40) 

and  from  equations  2-31  and  2-36 


1 


S.sm$  dO 


sin 


( ~iHj  jg,  iHs  f^T  2811 

\-ir‘ 


cosO 


(2  -  41) 


Solving  simultaneously  equations  2-37,  2-38,  2-39,  2-41  yields  equation  2-27.  As 
expected,  in  the  limit  that  the  skin  depth  is  larger  than  the  wall  thickness  the  low  frequency 
internal  magnetic  field  employed  in  equation  1-2  of  Appendix  A  is  obtained. 

Following  the  low  frequency  treatment  in  Appendix  A  (equation  A- 12)  the  aver¬ 
age  field  in  the  walls  is 


and  the  electric  field  at  the  wire  which  runs  in  th  ^  y-direction  (figure  7)  with  the  coordinates 
X  =  $  and  z  =  Tj  is 


3Ay  ^  /  c\ _ _ ^ 

dt  \  2/  co8(j9A)  —  ^  sin(/?A)  b  +  c 


(2  -  43) 
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2.^  Intermediate  Frequency  Result 


The  wire  current  follows  from  substituting  E*'  and  E*  from  equations  2-25  and  2- 
43  in  equation  2-6.  Because  the  electric  field  E”  created  by  the  external  electromagnetic 
pulse  is  uniform  along  the  wire  it  drives  a  uniform  current  down  the  wire  and,  consequently, 
only  the  m  =  0  mode  in  equation  2-25  is  retained.  The  wire  current  is,  therefore. 


io/TioH  [t]  2^cos(/9A)  —  |/?Ssin()9A)  b  +  c 

(1  -  i)  .8u;';^r^in»(^)sin^(ap) 

2?rwoVo  (;)*  +  (?)* 


(2  -  44) 


Clearly,  the  real  and  imaginary  parts  of  the  denominator  can  be  considered  to  be  the 
resistive  and  inductive  parts  of  the  wire  impedance  within  the  box. 

The  only  significant  approximation  made  in  the  derivation  has  been  to  approx¬ 
imate  the  internal  field  by  the  field  of  a  sphere.  Appendix  A  indicates  that  this  is  a  very 
good  approximation  in  the  small  number  of  cases  which  are  treated.  Perhaps  the  most 
convincing  argument  that  the  approximation  is  always  valid  lies  in  figure  8  which  suggests 
that  the  internal  field  is  insensitive  to  the  enclosure  geometry. 
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SECTION  3 


HIGH  FREQUENCY-PARALLEL  PLATE  ENCLOSURE 


3.1  Introduction 

At  ver}’  high  frequencies  not  only  is  the  skin  depth  smaller  than  the  wall  thickness 
but  the  electromagnetic  modes  of  the  enclosure  are  excited.  Consequently,  the  internal 
magnetic  field  can  no  longer  be  expected  to  be  constant.  In  sections  4  and  5  the  problem 
is  treated  precisely  but  the  formal  solution  is  complex  and  somewhat  inconvenient  to  use. 
In  this  section  the  physics  of  the  high  frequency  regime  is  examined  by  calculating  the 
current  excited  on  a  wire  within  the  enclosure  formed  by  two  infinite  parallel  conductors 
(figure  10).  The  simplified  geometry  allows  an  exact  solution  to  be  derived. 


3.2  Parallel  Plate  Enclosure 


Figure  10  illustrates  an  enclosure  formed  by  two  parallel  conducting  plates.  On 
the  external  surface  of  both  plates  is  an  oscillating  electric  and  magnetic  field 


Hx  =  Hie-‘“'‘ 

Ey  =  (3-1) 

as  illustrated.  These  fields  might  be  produced,  for  instance,  if  each  plate  were  exposed 
to  a  plane  electromagnetic  wave.  The  fields  Hi  and  Ei  are  the  sum  of  the  incident  and 
reflected  wave.  Again  the  wire  current  must  satisfy  equation  2-6  and  it  is  only  necessary 
to  calculate  E'^  and  E*. 
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wire 


Figure  10.  Parallel  plate  enclosure. 
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Calculation  of  E* 

If  the  wire  is  absent  then  the  electric  and  magnetic  fields  within  the  enclosure 
satisfy  the  free  space  wave  equations.  Because  the  cavity  excitation  is  identical  on  both 
plates  the  internal  magnetic  field  is  symmetrical  about  z  =  0.  Hence,  the  internal  field  is 


H'  =  H4C08(kz)e 


—  iult 


(3-2) 


where  k  =  ujv.  The  ini;  rual  electric  field  can  be  found  by  Ampere’s  law 


„  1  dH 

Ey - —  .  - 

IWCo  oz 


(3-3) 


to  be 


— iu»t 


Ev  =  : — H^sinfkzle 

»6oV 


Within  the  conductor  the  two  independent  solutions  of  the  wave  equation  are 

Hje‘^* 


(3-4) 


and 


(3-5) 


and  the  electric  field  of  these  modes 


and 


(3-6) 
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can  be  calculated  from  Faraday’s  law. 

Equating  the  tangential  components  of  E  and  H  at  z=a  and  the  tangential  com¬ 
ponents  of  H  at  z=b  yield  three  algebraic  equations  which  can  be  solved  for  the  unknowns 
H],  H3  and  H4.  Solving  yields  the  internal  magnetic  field  in  terms  of  the  external 


_ Hi _ 

co8(ka)  cos(/?A) - sin(ka)  sin(;3A) 

cov/J 


(3-7) 


From  equation  3-4  the  internal  electric  field  is 


E,  = - 

cos(ka)  cos(/?A)  —  sin(ka)  sin(/?A) 

As  expected,  equation  3-7  reduces  to  equation  2-28  in  the  limit  that  the  frequency 
is  so  low  that  kaC  1.  Comparing  these  expressions  indicates  that  the  high  frequency  modal 
structure  of  the  enclosure  enters  only  through  the  sin  and  cosine  of  the  argument  ka. 


3.3  Calculation  of  Wire  Driver 


The  coordinate  system  in  figure  11  is  employed  to  calculate  E'*  the  field  produced 
by  the  wire  if  the  external  driving  wave  is  absent.  The  parallel  walls  are  again  approximated 
as  perfect  conductors  which  is  an  excellent  assumption  at  high  frequency  when  the  walls 
are  good  shields.  The  walls  are  located  at  z  =  0  and  z  =  2a  and  the  wire  runs  parallel  to 
the  y  axis  and  is  located  at  x  =  0  and  z  =  $. 

The  case  is  treated  in  which  the  wire  current  and  electric  field  are  independent 
of  the  y-  coordinate.  The  y-  component  of  the  electric  field  obeys  the  wave  equation 


a*E 

dx^ 


+ 


d^E  u} 


dz^ 


+  -tE  =  -iw^oJ 


(3-9) 


subject  to  the  boundary  condition  that  it  vanishes  at  the  walls.  This  equation  will  be 
tran.sformed  by  the  fourier  sine  series 
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Figure  11.  Parallel  plate  enclosure. 
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En(x)  =  -/^  E(x,z)sin— da 

E(x,*)  =  EEnW*!”^  (3-10) 

nal 

in  the  z  variable  and  fourier  transformed  in  the  x  variable  by 


Mk)  =  E.(x).-'^*c 

E.(x)  =  /E„(ky‘»dk  (3-11) 

Then  the  transformed  electric  field  is 


„  „ ,  iwMoJn(k) 

En(kJ  =  — j[w 


k»  +  (~)‘  -  (f  )* 


(3  -  12) 


where  Jn(k)  is  the  transformed  wire  current  density  J.  The  convolution  theorem  indicates 


E.(x)  =  ^  /  dxJ„(x)f(x  -  x')  (3  -  13) 

where 


dke‘^ 


WfiQ 

(k  -  a){k  +  a) 


(3  -  14) 


and 


(3  -  15) 


f 
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It  is  assumed  that  the  field  damps  very  slightly  due  to  an  unspecified  mechanism.  Then  u 
must  have  a  small,  negative  imaginary  part  since  the  temporal  dependance  is  exp  (-i  uf  t). 
The  poles  are  at  a  and  -a.  The  square  root  appearing  in  the  definition  of  a  is  chosen  to 
be  that  branch  of  the  function  for  which  if 

u  =  |u|e'^  —  IT  <  9  <  ir 


then 


(3  -  16) 


Figure  12  shows  u/v  with  its  small  imaginary  part.  Also  indicated  is  the  value 
of  a  for  several  values  of  the  integer  n.  The  corresponding  values  of  -a  are  shown.  The 
inverse  transform  is  done  by  closing  the  contour  in  the  upper  k  plane  when  x  is  greater 
than  zero.  Yielding 


Finally,  the  inverse  sine  series  is  inverted  to  yield  the  desired  value  of  S'* 


E*(x,z) 


ujfiol  ^  .  ,n7r$  e 


VO’ -O' 


.  ,n7rz 
«n(— ) 


(3-17) 


The  field  at  the  wire  is  found  by  evaluating  this  expression  at  x  =  0  and  z  =  It  is 
assumed  that  the  wire  current  density  is  J  =  I  5(z  -- 


3.4  Wire  Current  in  Parallel  Plate  Enclosure 


The  wire  current  is  found  to  be 
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31  » 


-i/iovHi  sin(ki/) 


2e‘^^ 


1  = 


COB  ka  +  sin  ka 

JL 


nail  ^  ^ 


(3  -  18) 


by  substituting  equations  3-8  and  3-17  in  equation  2-6.  By  definition  0A  in  the  high 
frequency  regime  is  large  so  the  large  argument  limits  of  the  sin  and  cosine  functions  have 
been  employed.  Furthermore,  equation  3-17  has  been  adjusted  to  conform  to  the  figure  10 
coordinate  system  by  evaluating  the  wire  location  at  i/  =  (  +  a. 

At  first  glance  it  might  be  thought  that  the  wire  current  would  be  very  large  if 
the  incident  wave  frequency  equals  the  frequency  of  one  of  the  modes  of  the  cavity.  This 
analysis  indicates  that  the  wire  current  is  large  but  not  huge  at  the  modal  frequencies. 
The  electric  field  induced  within  the  cavity  is  largest  when  the  expression  (equation  3-18) 


cos(ka)  -I-  sin(ka)  (3  -  19) 

is  smallest.  For  practical  values  of  the  conductivity  the  minimum  value  is  one  which  occurs 
when  sin  ka  =  0.  As  expected  the  minimum  occurs  whenever  the  driving  frequency  equals 
the  frequency  of  one  of  the  internal  modes  of  the  cavity  (if  its  wall  conductivity  were 
infinite).  But  is  a  complex  number.  Hence,  there  is  no  real  value  of  k  =s  cj/v  for  which 
equation  3-19  vanishes.  Therefore  the  internal  field  is  finite.  Evidently,  if  the  walls  are 
made  of  metal  then  the  conductivity  always  introduces  enough  damping  to  prevent  the 
internal  field  from  becoming  very  large.  If  the  walls  were  made  of  a  very  poor  conductor, 
however,  then  equation  3-18  shows  that  the  internal  field  could  be  driven  to  large  values. 
The  upward  spikes  in  figure  5  illustrate  the  increased  current. 

The  square  root  in  equation  3-18  is  also  zero  whenever  the  driving  frequency 
equals  the  frequency  of  an  internal  mode.  This  causes  the  wire  current  to  be  zero  which 
is  not  surprising  for  at  for  these  frequencies  a  finite  wire  current  causes  an  infinite  electric 
field.  If  the  finite  wall  conductivity  were  taken  into  account  then  the  current  would  be 
small  but  not  zero.  Consequently,  equation  3-18  indicates  that  the  wire  current  at  these 
frequencies  is  zero  which  is  an  underestimate.  The  downward  spikes  in  figure  5  illustrates 
this  behavior. 

The  internal  field  is  driven  maximally  whenever  ka  =  n^r  or 
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(3  -  20) 

whereas  the  square  root  is  zero  if 

ItiK  1 

—  =  -mjT  .  (3  -  21) 

The  difference  is  attributable  to  the  fact  that  the  external  field  is  assumed  to  be  equal  on 
each  conducting  plate.  Hence,  it  can  excite  only  modes  of  E*  which  are  anti-symmetric 
in  z.  The  wire,  however,  excites  both  symmetric  and  anti-synunetric  modes  depending  on 
its  location.  It  should  be  pointed  out  that  the  peak  current  at  each  upward  spike  is  also 
underestimated  because  whenever  equation  3-20  is  satisfied  equation  3-21  is  satisfied  also 
for  m  =  2n. 


b>a 

—  =  nx 

V 
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SECTION  4 


HIGH  FREQUENCY  REGIME  -  APPROXIMATE  TREATMENT 


In  the  low  and  intermediate  frequency  regime  the  electromagnetic  field  within 
the  enclosure  is  approximately  uniform  and  easily  calculated  with  equation  2*27  and  2-43 
of  section  2.  At  high  frequencies,  however,  the  internal  field  depends  upon  the  frequency 
and  angle  of  incidence  of  the  incident  electromagnetic  wave  which  greatly  complicates  the 
problem  of  predicting  the  wire  current.  Several  approximations  are  made  in  this  section 
in  order  to  derive  an  easily  applied  formula  to  estimate  the  wire  current.  In  section  5  a 
more  exact  treatment  is  derived. 

4.1  Inductance  Matrix 

In  this  section  electromagnetic  transmission  through  a  conducting  plane  is  ex¬ 
amined  in  some  detail.  Figure  13  illustrates  a  conductor  of  thickness  A  which  is  irradiated 
normally  from  the  right  with  a  linearly  polarized,  plane  electromagnetic  wave  of  frequency 
u.  The  total  field  at  z  =  A  is  the  sum  of  the  incident  and  reflected  waves  and  is  denoted 
by  E‘  and  Similarly,  the  total  field  at  the  z  =  0  plane  is  and  H*. 

Within  the  conducting  medium  two  independent  solutions  to  the  diffusion  equa¬ 
tion  exist.  The  total  internal  magnetic  fields  is 

H(z,  t)  =  (4  -  1) 

The  total  electric  field  is  derived  from  Faradays  law  to  be 

E{z,t)  =  — Hje‘^‘  -  — Hse-‘^*  (4  -  2) 

<7  O 

The  wave  H*  travels  to  the  right  while  H’  travels  to  the  left. 

The  electric  and  magnetic  fields  are  continuous  across  the  surface.  Equating  the 
internal  and  external  fields  at  the  two  surfaces  leads  to  four  algebraic  equations  which 
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contain  the  six  quantities  through  and  and  E^.  It  is  useful  to  consider  and 
to  be  known  and  solve  for  the  remaining  quantities  in  terms  of  them.  Then 


E,  =  Z,Hi-Z,H4 


E4  =  ZtHi-Z,H4 


(4-3) 


where 


Z.  = 


Zt 


P  cos/^A 
a  sin/?A 

/? 

ffsin/JA 


(4-4) 


are  the  self  impedance  and  transfer  impedance.  This  can  be  written  in  matrix  form  eis 


E,  ■ 

■  Z. 

-Zt 

Hi  ■ 

.  E4. 

. 

-Z. 

.^4  . 

(4-5) 


In  terms  of  the  fields  on  the  z  =  A  surface,  and  E^  the  magnetic  field  phasors 
in  the  metal  are 


nr 

Hj  =  ,  (4-6) 

The  electromagnetics  are  explained  in  figure  14  assuming  that  the  driving  wave  is 
incident  at  the  z  =  A  surface.  On  the  z  =  0  surface  the  ratio  of  electric  and  magnetic  fields 
within  the  conductor  is  constrained  to  be  equal  to  the  impedance  Z,  of  the  electromagnetic 
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Figure  14. 


Qualitative  behavior  of  the  electric  and  magnetic  fields  in  figure 
13  in  the  high  frequency  regime. 
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field  just  oiitside  the  metal.  To  fulfill  this  requirement  the  II*  wave  is  sent  backward  from 
the  z  =  0  surface  into  the  conductor.  Usually,  the  two  waves  in  the  metal  have  nearly  equal 
and  opposite  values  of  H  on  the  z  =  0  surface  in  order  to  reduce  the  total  H  field.  The 
reduction  is  necessary  because  the  conductor  supports  waves  with  much  larger  ratios  of  H/E 
than  are  usually  found  in  vacuum.  The  figure  indicates  that  the  backward  wave  magnitude 
is  very  small  in  the  high  frequency,  small  skin  depth  regime  so  that  approximately 


H3e-‘^^  =  Hi  . 


Continuity  of  the  electric  and  magnetic  field  at  z  =  0  yield 


(4  ™  7) 


H4  =  Hj  +  Ha 


a 


(4-8) 


Substituting  equation  4-7  in  4-8  and  solving  for  the  internal  electric  field  yields 

+  — H4 
a  a 

which  is  the  same  as  the  more  precise  equation  4-3  in  the  high  frequency  limit  where 
|/?A|  »  1. 

But 


K 


(4-9) 


so  that  equation  4-3  yields 


E4  = 


ZtH, 


1  +  Z,/Ze 


(4  -  10) 


In  many  cases  Zg/Z*  ■<!.  In  the  parallel  plate  czise,  for  instance,  equation  3-4  in  section  3 
shows  that 
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Z«  =  i/iov  sin  ka 


(4  -  U) 


ajid  the  inequality  is  true  provided 


sinka 


p  cosfidk. 
i/iova  sin/9A 


% 

(4  -  12) 


where  k  =  w/v. 

But  sin  ka  =  0  is  the  resonance  condition  for  the  inodes  in  a  parallel  plate 
enclosure  (with  infinitely  conducting  v/alls).  So  within  a  very  narrow  frequency  band 
about  each  internal  mode  equation  4-12  is  not  fulfilled.  Equation  4-12  indicates  the  width 
of  the  band  (Acv)  is 


Atv  _  6 
cvo  a 

where  wq  is  the  resonance  frequency  given  by  sin( - )  =  0  and  6  is  the  skin  depth  in  the 

walls.  At  all  other  frequencies  the  self  impedance  term  in  equation  4-3  can  be  neglected 
and  the  internal  electric  field  is  very  well  approximated  by 


E4  =  ZtHi  .  (4  —  13) 

Within  the  equation  4-12  frequency  band,  however,  the  impedance  of  the  enclo¬ 
sure  and  the  impedance  of  the  conductor  are  very  well  matched  so  that  the  magnetic  field 
H4  within  the  enclosure  is  comparatively  large.  This  very  narrow  resonance  was  pointed 
out  in  the  section  3  discussion  of  the  parallel  plate  enclosure.  The  method  derived  in 
the  remainder  of  this  section  employs  the  equation  4-13  simplification  and  is,  therefore, 
not  valid  in  the  narrow  frequency  band  about  each  mode  of  the  enclosure.  In  section  5  a 
method  is  derived  to  treat  the  narrow  band  of  frequencies. 
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4.2  Calculation  of  £* 


Figure  15  illustrates  the  box  with  the  wire  running  in  the  y  direction.  The  Ey 
component  of  the  electric  field  obeys  the  wave  equation 

V*Ey  +  k’Ey  =0  (4  -  14) 

Section  A  points  out  that  in  the  absence  of  the  wire  the  solution  to  this  equation  can  be 
written  in  terms  of  surface  integrals  over  the  walls.  It  is 


where  G 


E, 


i.c  wall! 


-  i;;  / 


(4-15) 


.  ,£tx.  /tnjry.  .  ,n7rz.  ,  ,t7rx\  ,Tniry' .  ,  ,nxz'. 

»  ,  sm(— )co4(-^)sin(— )3m(— )c.,s(— )s.n(-^)  1 

Zw  g,bc  r  .  »  1  UJ 


£«nsl 

m=0 


(^)’  +  (r)’-(7r 

a  D  c 


k* 


(4  -  16) 


is  the  Green’s  function  which  satisfies  Dirichet  boundary  conditions  on  the  x  =  constant 
and  z  =  constant  wails  and  Neumann  boundary  conditions  on  the  y  =  constant  walls. 


The  integrals  contain  Ey  on  the  x  ~  constant  and  z  =  constant  walls.  But 
equation  13  of  the  last  section  shows  that  these  quantities  are  easily  calculated.  On  the  x 
=  constant  walls  they  are 


Ey(x  =  0,y,z)  =  +ZtH, 

Ey(x  =  a,y,z)  =  -ZtH,  (4-17) 
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and  on  the  z  =  constant  walls  they  are 


Ey(x,y,z  =  0)  =  -Z,Hx 

Ey(x,y,z  =  c)  =  +ZtHx  (4-18) 

where  the  magnetic  intensities  H  are  on  the  external  surfaces  of  the  box.  It  is  intended  that 
the  external  magnetic  fields  will  be  calculated  with  the  code  HFC  which  is  explained  in 
section  6.  HFC  assumes  the  box  has  infinite  conductivity  and  calculates  H  on  the  external 
surface.  This  is  an  excellent  approximation.  Other  approximations  for  H  on  the  external 
surfaces  which  avoid  the  use  of  a  code  are  mentioned  below. 

On  the  y  =  constant  surfaces  the  integral  requires  the  normal  derivative  of 
dEy/dy  at  the  surface.  A  small  distance  from  the  y  =  constant  walls  (inside  the  box) 
the  divergence  of  the  electric  field  is  zero  so  that 

aEy  ^  ^ 

dy  dx  d% 

On  the  y=b  surface,  for  instance,  the  right  hand  member  is  derived  from  equation  4-13  to 
be 


aEx  _  „  aH, 
ax  ax 

(4- 

where  Hx  and  H.  are  the  magnetic  intensity  on  the  external  surface  of  the  y=b  wall.  They 
could  be  supplied  by  the  code  HFC.  Consequently,  equation  4-15  yields  the  internal  field  Ey 
in  terms  of  a  quadrature  over  the  external  magnetic  field.  Of  course,  it  has  been  assumed 
that  the  one  dimensional  transfer  impedance  in  the  previous  section  applies.  This  is  well 
justified  since  the  walls  are  much  thinner  than  cither  the  wavelength  or  the  lengths  of  the 
box’s  sides.  Within  a  very  narrow  frequency  band  around  each  of  the  box’s  modes  this 
method  fails.  Section  5  explains  how  to  calculate  the  fields  within  these  bands. 


dE.  _  ^ 

dz  dz 
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In  the  most  general  case  the  external  field  is  incident  on  the  box  at  arbitrary 
angle  and  may  have  any  spatial  dependence  satisfying  the  wave  equation.  HFC  can  treat 
this  case  and  will  furnish  the  magnetic  intensity  H  on  the  external  surface  of  the  box. 
Equations  4-17^  4-18  and  4-19  then  furnish  Ey  and  dEy/dy  from  H.  Finally,  equation  4-15 
yields  Ey  inside  the  box  which  is  needed  in  the  next  section  to  calculate  the  wire  current. 

The  integral  in  equation  4-15  is  reasonably  straight-  forward.  In  some  cases  H 
approaches  infinity  at  the  edge  of  the  box.  But  these  are  integrable  singularities.  In  order 
to  illustrate  the  method  a  simple  example  is  worked  out  in  the  next  subsection. 


4.3  Simple  Example 

In  order  to  illustrate  the  general  method  in  the  previous  subsection  consider  the 
case  of  a  plane  wave  normally  incident  on  the  z  =  c  wall  of  the  box  in  figure  15.  If  the 
wavelength  is  much  less  than  the  box  dimensions  then  the  geometric  optics  theory  of  light 
indicates  that  total  magnetic  field  on  the  z  =  c  face  is  twice  that  of  the  incident  field  and 
is  constant  over  the  surface.  The  surface  magnetic  field  is  smaller  on  the  other  faces.  In 
this  example  it  is  assumed  that  the  surface  magnetic  field  is  zero  everywhere  except  the 
exposed  face  and  is  well  behaved  at  the  edges  of  the  box. 

Equation  4-18  shows  that  the  internal  electric  field  on  the  z  =  c  wall  is 

Ey  =  2ZtH;  (4  -  21) 

where  is  the  field  of  the  incident  wave  generated  by  a  source  of  electromagnetic  inter¬ 
ference,  perhaps.  The  electric  fields  tangent  to  the  other  walls  will  be  much  less  than  this. 
Equation  4-3  shows  that  these  electric  fields  can  be  computed  from  Hi  and  H4  the  external 
and  internal  magnetic  fields  on  each  v/all.  But,  by  assumption  Hi  is  zero  except  on  the 
z=c  face.  Consequently,  the  internal  electric  field  on  the  other  walls  is  given  by 


E4  =  -Z,H4  (4  -  22) 

The  internal  magnetic  field  H4  is  roughly  the  same  on  all  walls  of  the  box.  Hence,  the 
magnitude  of  the  internal  field  E4  is  about  the  same  on  all  the  walls.  But  in  section  4.1 
this  term  was  shown  to  make  a  negligible  correction  to  the  internal  electric  field  on  the  z  == 
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c  wall.  Hence,  the  internal  held  is  much  larger  on  the  exposed  wail  (z=c)  than  the  others. 
From  equation  4-19  it  can  be  inferred  that  dEy/dy  is  small  also. 

Consequently,  only  the  integral  over  the  z  =  c  face  survives  in  equation  4-15  and 


E,(x,z)  = 


8Z.HI 

7r*c* 


E 


<,n=l 


ncosn7r(l  —  cos^x) 

<!(;)’  + 

a  c 


.  ^nx  . 

sm - sin 

a 


Ey  is  independent  of  y  for  this  case. 


(4  -  23) 


4.4  Wire  Current  E*' 


At  each  point  on  the  wire  the  relationship 


IRo(l-i)^  =  E'  +  E*  (4-24) 

must  be  satisfied.  It  is  convenient  to  transform  each  side  of  this  equation  with  a  cosine 
transform 


yielding 


fm  = 


fb 

i 


miry 

cos  — T — dy 


^(y) 


00 

fm  cos 

m=:0 


miry 

~ 


I,„R„(1  -  i)^  =  e;  +  E 


where  the  subscript  m  indicates  the  ni*'*  fourier  coefficient. 


(4-25) 


(4  -  26) 
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Simple  Case 

First  the  wire  current  will  be  derived  for  the  simple  case  treated  in  the  previous 
section.  Because  the  plane  wave  is  normally  incident  on  the  box,  Ey  (equation  4-23)  is 
independent  of  the  y  -coordinate  and  only  the  m  =  0  coefficient  is  non-zero.  It  is  trivially 
given  by  equation  4-23  to  be 


SZtHj. 

f,n=l 


n(cosn7r)(l  —  cos/tt) 


.  .  ,/mc.  .  .nxz. 
)sin(-— )sm(— ) 
a  c 


(4  -  27) 


In  section  A. 2  the  electric  field  E*  produced  by  the  wire  was  calculated.  From 
equation  2-25  it  is  clear  that  the  m  =  0  component  is 


e:=o  =  -7 


.  5  .  I 

A,  ,T  oo  Sin  - sm  - 

-dwlmso  r-'  a  r 


k*  -k* 


Substituting  equations  4-27  and  4-28  in  4-26  yields  the  wire  current 


(4  -  28) 


I  = 


-8ZtHi  ^  f  n(cosn7r)(l  -  cos£7r)  I  . 

■^£.1  <i(^)+Tf)~kTr‘ 


.  n-nri 

sin - sm - 

a  c 


..Wo  4ia; 


^  ain^^sin^gf 


2^0  aceov*  k*  -  k* 


(4  -  29) 


56 


General  Case 


For  an  arbitrary  incident  wave  the  code  KFC  (section  6)  furnishes  the  magnetic 

intensity  H  on  the  external  surface.  Equations  4-17,  4-18  and  4-19  then  yield  the  internal 

electric  field  on  all  surfaces  which  are  used  in  4-15  to  yield  the  internal  electric  field  at  the 

niTry 

wire  E*'.  The  Green’s  function  in  equation  4-15  contains  the  factor  cos  — ; — .  Hence,  this 

b 

formalism  will  always  furnish  E*  as  a  cosine  transform  of  the  form 


E*  = 


OO 


E  Em  cos 


m^y 

T" 


Of  course,  E*  will  contain  the  integrals  in  equation  4-15  but  these  are  straightforward  to 
perform  numerically  in  the  general  case. 

From  2-25  the  m***  fourier  coefficient  of  E*’  is 


E! 


—41  “ 

41.  ^ 


-(-y 


iwacco  k*  -  k» 


.  2  .  2 

sin - sin  - 

a  c 


Equation  4-26  then  yields  the  m^**  Fourier  coefficient  of  the  wire  current 


I.  = 


—  E 


^(1  •" 


._2  _._2  OTTl? 

k'-k*  a  C 


Finally,  the  current  1  in  the  wire  at  any  position  is  given  by  summing  the  Fourier  compo¬ 
nents 


1  = 


OO 

H  In.  cos 

m=0 


mTry 

“b~ 


The  calculation  is  straightforward,  but  tedious.  It  is  expected  to  be  valid  except 
in  the  narrow  frequency  band  about  each  internal  mode  of  the  box.  The  behavior  of  the 
fields  at  the  edges  of  the  box  should  be  considered  more  thoroughly. 
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SECTION  5 


HIGH  FREQUENCY  REGIME  -  COMPLETE  TREATMENT 


It  has  been  pointed  out  that  the  section  4  method  does  not  apply  in  a  very  narrow 
frequency  band  centered  about  each  mode  of  the  enclosure.  Within  the  band  the  internal 
magnetic  field  (equation  4-3)  is  driven  to  comparatively  large  values  and  the  internal 
electric  field  is  not  accurately  approximated  by  ZtHi.  The  internal  electric  and  magnetic 
fields  must  be  calculated  self  consistently.  A  method  of  calculating  the  wire  current  in  this 
situation  is  outlined  in  this  section. 

Consider  the  component  of  the  internal  magnetic  intensity  which  in  the  ab¬ 
sence  of  the  wire  satisfies  the  homogeneous  Helmholtz  equation 

V*H.  -h  =0  (5-1) 

In  terms  of  the  magnetic  intensity  at  the  inside  surface  of  the  walls  section  A 
points  out 


H,  =  /  dS' .  H,  V'G  (5  -  2) 

47r  J 

where  the  surface  integrals  are  taken  over  the  x  =  constant  and  y  =  constant  walls.  Because 
He  is  normal  to  the  z  =  constant  walls  it  is  very  small  on  these  surfaces  and  is  neglected. 
The  Green’s  function  satisfies  Dirichlet  boundary  conditions 


^  / abc . 

G  =  47r(— ) 


,  tirx  .  mjry  .  nxz  .  tirx.'  .  m7ry'  .  n^rz' 
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. a _ b _ c  a _ b _ S_ 


<,in,n=l 
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(5-3) 


into 


The  electric  field  tangent  to  these  four  interior  walls  can  be  fourier  decomposed 
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Ey(0,y,z)  =  X) 


2y0 


m,n=l 


Ey(a,y,2)  = 


in,n— 1 


^  .  mTry  .  nirz 
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'“be 


Ex  (x,  0,2)  = 
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2-  e^sin  — sm  — 


Ex(x,b,z)  = 
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a  c 


(5-4) 


Equation  4*3  indicates  that  the  the  x  =  0  wall 


H.  =  ^Ex  +  2Hie‘'»^  (5-5) 

where  Hi  is  the  external  held  on  the  surface  of  the  box  which  is  computed  by  the  code 
HFC  (section  6).  On  the  x  =  a,  y  =  0,b  walls  analogous  equations  furnish  the  interior 
surface  electric  fields. 

In  addition,  on  the  y  =  b  wall  Faraday’s  law  shows 


-47r  ^H, 
iwco  dy 


Substituting,  equations  5-4,  5-5,  5-6  into  5-3  yields 


(5-6) 
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l-K  .  a  be  _ 
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(5-7) 


Here  the  Fourier  eomponents  on  the  x=a  wall  of  the  field  Hi  is  denoted  by  H^"(x  =  a). 
Other  components  of  Hi  are  denoted  by  a  similar  notation. 

This  equation  relates  e™,e™,e^,e^  and  the  H^  components  given  by  the  code 
HFC  on  the  x  =  0,a  and  y  =  0,b  walls.  Equations  5-4  have  4S*  unknown  Fourier  coefficients. 
Equation  5-7  is  S*  equations  relating  these  unknowns.  3S*  more  equations  are  found  by 
applying  the  same  method  to  calculate  H,  on  the  y  =  0  and  x  =  0,a  planes.  Solving  these 
simultaneous  algebraic  equations  numerically  furnishes  all  the  coefficients  in  equations  5-4. 
Some  care  must  be  exercised  when  evaluating  the  series  expressions  near  the  edges  since 
they  may  converge  slowly. 

Applying  the  same  technique  to  the  Hx  and  Hy  counterparts  of  equation  5-1 
allows  a  self  consistent  calculation  for  the  internal  electric  field  tangential  to  all  the  internal 
surfaces  of  the  box.  From  these  E*  can  be  computed.  Section  2  yields  E'*  which  can  then 
be  employed  with 


IRo(l-i)~  =  E*4-E'^ 

to  furnish  the  wire  current  as  in  sections  2,  3  and  4. 
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SECTION  6 


HIGH  FREQUENCY  CODE 


A  rectangular  box  illuminated  by  electromagnetic  radiation  can  be  described  in 
terms  of  the  magnetic  fields  or  alternatively  in  terms  of  the  surface  current  density  K.  A 
careful  analysis,  yields  the  equation:’ 


K(f)  -  2y'n  X  K  X  G(r,f')dS'  =  2  n  x  H«*(f)  (6  -  1) 

where  K  is  the  surface  current,  H***  is  the  external  magnetic  field  and  G  is  the  outgoing 
Green  function 


G(f,?) 


_ 1 _ 

47r|r  —  ?|exp(ik|r  —  Fj) 


(8-2) 


corresponding  to  frequency  u  =  kV.  The  integral  over  S'  is  over  all  the  faces  of  the  rectan¬ 
gular  enclosure.  Note  that  only  the  tangential  component  of  the  incident  field  contributes 
to  the  formation  of  current  on  the  surface.  The  enclosure  walls  are  infinitely  conductive. 

In  order  to  solve  the  integral  equation  (6-1)  for  K  we  have  written  a  computer 
code  called  High  Frequency  Code  (HFC).  This  code  solves  the  integral  equation  above 
numerically  using  finite  difference  approximation  to  the  double  integral.  In  principle,  the 
code  permits  as  fine  a  resolution  of  the  surface  as  necessary  because  of  the  utilization  of 
virtual  memory  in  the  computer.  In  practice,  the  resolution  is  limited  by  both  the  accu¬ 
racy  that  can  be  obtained  inverting  a  large  matrix  and  by  the  amount  of  computer  time 
available.  Thus  accuracy  limits  the  resolution  to  36  elements  per  face  in  single  precision 
arithmetic  and  to  about  200  elements  per  face  in  double  precision  arithmetic.  The  time 
required  for  computation  increases  as  the  third  power  of  the  total  number  of  surface  ele¬ 
ments  used.  For  a  4  MIPs  machine,  100  elements  per  face  requires  about  1  hour  of  CPU 
time  which  costs  about  $20  on  an  ELXSI  minicomputer. 

In  general  HFC  can  treat  any  angle  of  incidence  and  any  frequency.  Round  off 
error,  however,  may  prevent  HFC  from  treating  frequency  greater  than  about  ten  gigahertz 
for  boxes  about  one  meter  on  each  side. 
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For  the  present  problem  we  have  considered  only  a  rectangular  box.  Other  box 
geometries  could  be  also  computed  but  the  algebra  becomes  somewhat  more  involved.  For 
the  six  faces  of  the  box,  there  are  12  independent  components  of  the  current  together 
with  the  boundary  conditions  at  the  edges  of  the  box  which  arise  from  continuity.  Each 
component  is  complex  and  the  physical  current  is  the  real  part  of  the  product  K  exp(iu;t). 

The  resultant  matrix  equation  can  be  cast  in  the  form: 


AK  =  D 


where  A  is  the  LHS  matrix  (kernel)  and  D  is  the  driving  term.  In  order  to  save  memory, 
we  had  hoped  that  the  kernel  has  some  evident  symmetries  that  permit  more  efficient 
computation.  However,  as  it  turns  out,  the  matrix  A  is  band  structured,  but  not  evidently 
symmetric  (or  antisymmetric).  Thus  the  computation  of  the  surface  currents  requires  large 
memory  available  unless  we  invoke  the  physical  symmetries,  which  though  not  apparent 
must  be  present  in  the  solution.  These  include  rotational  and  transinational  symmetries 
which  are  clearly  embedded  in  the  original  Maxwells  equations.  Our  calculations  indicate 
that  these  are  indeed  satisfied  by  our  numerical  solutions. 

To  test  the  validity  of  the  code  we  have  computed  the  surface  currents  on  a  long 
rectangular  box.  For  this  situation  and  an  incidence  perpendicular  to  the  long  axis  of  the 
box,  the  solution  has  been  computed  previously  e.g.  by  Mei  and  van  Bladel.®  Their  results 
together  with  our  computations  are  shown  in  figure  16.  In  view  of  our  finite  resolution  in 
the  finite  difference  scheme  the  agreement  with  their  results  is  excellent.  The  results  are 
somewhat  less  accurate  for  the  E  wave,  i.e.  one  where  the  magnetic  field  is  perpendicular 
to  the  long  axis  of  the  cylinder.  This  is  because  the  surface  current,  and  the  induced 
magnetic  field  is  now  discontinuous  and  divergent  the  edges  perpc:ndicular  to  the  magnetic 
field.  This  can  be  seen  in  figure  17  where  we  show  results  from  reference  8  for  this  case.  We 
have  attempted  to  compute  this  problem  with  our  code  with  some  success,  but  to  obtain 
accurate  results  requires  more  time  than  we  have  presently  available. 


6.1  CODE  DESCRIPTION 


The  basic  tenet  that  has  been  I'sed  in  designing  the  HFC  code  was  to  maincam 
minimum  complexity,  in  particular  since  code  development  was  not  our  primary  goal  for 
this  effort.  For  this  reason  the  code  contains  minimum  “bells  and  whistles”  coding.  The 
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Figure  16.  Surface  current  on  a  long  square  box  in  an  H  wave.  The  triangles 
correspond  to  ka=l,  circles  to  ka=.3.  The  solid  lines  are  square 
cylider  calculations  of  Reference  8. 


Surface  current  density  for  a  square  cylinder  in  an  £  wave.  (Taken 
from  Refererence  8.)  Note  the  singular  behavior  near  the  edges. 


input  consists  of  the  dimensions  of  the  enclosure,  frequency  of  incident  radiation  and  the 
size  of  the  driver  field.  In  addition,  the  input  requires  the  number  of  vertices  in  each 
Independent  direction  x,y,z. 

The  code  consists  of  the  main  control  program,  which  defines  the  loop  structure 
necessary  to  compute  the  surface  integrals,  the  calls  to  the  algebraic  equation  solver  rou¬ 
tines  and  the  output  instructions.  The  code  uses  the  LINPACK  GESL  routine  to  solve  the 
algebraic  equations  for  K  in  complex  arithmetic.  A  GREEN  subroutine  is  used  to  compute 
the  green  function  and  its  derivatives  needed  in  equation  6-1. 

The  running  time  for  the  code  has  been  obtained  from  a  number  of  numerical 
experiments.  We  conclude  that  the  CPU  time  required, 

CPUTIM  =  TZERO(nx  *  ny  -f-  nx  *  nz  -t-  ny  *  nz)®, 
where  TZERO  —  1.6“^  seconds  for  our  4  MIP  ELXSI  computer. 
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APPENDIX  A 


INTRODUCTION 


A.l  Summary  of  Results 


Electromagnetic  radiation  can  diffuse  through  a  thin  walled  metallic  box  and 
couple  to  internal  circuitry.  At  high  frequencies  for  which  the  skin  depth  is  much  smaller 
than  the  wall  thickness  a  metallic  enclosure  is  a  good  shield.  At  low  frequencies,  however, 
the  incident  magnetic  field  penetrates  the  walls  easily  so  that  the  magnetic  field  within  the 
box  approaches  the  magnitude  of  the  external  field  radiating  the  enclosure.  This  report 
discusses  the  diffusion  cf  a  low  frequency  electromagnetic  signal  through  a  thin  walled 
box  that  contains  a  wire  which  is  connected  between  opposite  walls  of  the  parallelpiped 
(figure  A-1).  A  quantitative  analytic  theory  is  derived  to  predict  the  current  on  the  wire. 
Easily  applied  for  mulas  are  found  which  furnish  the  wire  current  as  a  function  of  time  due 
to  any  external  driving  pulse.  Frequency  domain  solutions  for  the  wire  current  are  also 
shown.  In  order  to  establish  the  theory’s  credibility  its  predictions  are  compared  to  the 
wire  current  calculated  by  the  code  BOX4. 

The  theory  developed  below  indicates  that  the  current  I(amperes)  in  the  wire  is 

given  by 


I(w) 


iptowH*  wL(L  4-  2x) 

R  -  \(jjZ  ^  2(w  +  L) 


(A-1) 


where  u  is  the  frequency  of  the  external  radiation  incident  on  the  enclosure,  R  is  the  wire 
resistance,  £,  is  the  wire  inductance  and  the  box  dimensions  and  wire  location  are  denoted 
by  w,L,x  as  in  figure  A-2.  The  permiability  of  space  is  and  i  =  v'^^.  H*,  the  magnetic 
field  intensity  within  the  box  is  approximated  by 


H* 


1 


1 

iuijlnga* 

3 


(A-2) 
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Figure  A-1.  The  incident  electromagnetic  wave  diffuses  through  the  box  and 
induces  a  current  on  the  internal  wire. 
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Wall  Currant 


Figure  A-2.  Problem  geometry. 


where  a  is  the  conductivity  of  the  walls,  and  A  is  the  wall  thickness.  The  external  magnetic 
held  radiating  the  box  is  assumed  to  be  the  real  part  of 


H°e-‘"‘ 


Because  the  wavelength  of  the  low  frequency  radiation  is  much  larger  than  the  box  is 
independent  of  position  in  the  vicinity  of  the  enclosure.  The  parameter  a  is 


a  = 


1/3 


in  terms  of  the  volume  V  of  the  enclosure.  MKS  units  are  used  throughout  this  report. 

Figure  A-3  compares  the  wire  current  predicted  by  these  expressions  with  the 
current  computed  by  the  computer  code  BOX4  discussed  below.  The  agreement  is  within  a 
factor  of  four  over  five  orders  of  magnitude  of  frequency.  At  one  hertz  angular  frequency  the 
code  and  analytic  method  differ  by  six  percent  and  at  100000  hertz  the  analytic  predictions 
are  a  factor  of  four  too  small.  In  this  example  the  box  is  two  meters  on  each  side  and  its 
walls  are  one  millimeter  thick.  The  wall  conductivity  is  a  =  10®  (mho/meter).  The  wire 
resistance  is  0.0366  ohms  and  it  is  one-half  meter  from  the  y'=0  and  z'=0  walls.  The  wire 
is  one  centimeter  in  radius  and  the  external  magnetic  field  is  perpendicular  to  the  wire 
and  has  a  magnitude  of  one  ampere/meter. 

Table  A-1  indicates  that  equations  A-1  and  A-2  predict  the  wire  current  at 
different  positions  of  the  wire  at  one  hertz  angular  frequency.  Agreement  is  within  twenty 
percent  when  the  wire  is  varied  between  five  centimeters  from  the  wall  and  five  centimeters 
from  the  center  of  the  box  which  is  2  meters  on  each  side.  All  other  parameters  are  identical 
to  the  figure  A-3  case. 

Table  A-2  indicates  that  at  100000  hertz  the  analytic  method  correctly  evaluates 
the  dependence  of  wire  current  on  wire  position.  As  the  position  is  varied  the  analytic 
prediction  is  always  about  a  factor  of  three  smaller  than  the  code.  The  analytic  method 
always  underpredicts  the  current  by  the  same  factor  at  high  frequency  regardless  of  wire 
position.  As  in  figure  A-3  the  disagreement  at  high  frequency  is  due  to  the  approximations 
made  in  the  derivation  equations  A-1  and  A-2. 

More  extensive  comparisons  are  necessary  to  evaluate  the  accuracy  of  the  analytic 
formulas  over  a  large  range  of  box  dimensions,  wire  parameters  and  wall  conductivity 
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I  Wire  Current  (Amperes) 


Figure  A-3.  Wire  current. 
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Table  A>1.  Wire  current  as  a  function  of  wire  position. 
Wire  Position  Wire  Current  Wire  Current 


X  (cm) _ Theory  (Amp)  Code  (Amp) 


5 

3.3  X10-® 

3.5  xlO-® 

10 

3.1  xlO-® 

3.3  xlO-' 

50 

1.7  xlO-* 

1.6  xl0-« 

95 

1.7  xl0-« 

1.4  xlO-® 

Table  A-2. 

Wire  current 

as  a  function  of  wire  position 

Wire  Position 

Wire  Current 

Wire  Current 

Ratio 

X  (cm) 

Theory  (Amp) 

Code  (Amp) 

Code  -r  Theory 

5 

0.033 

0.10 

3.0 

10 

0.027 

0.081 

3.0 

50 

0.011 

0.043 

3.9 

95 

0.0011 

0.0029 

2.6 

and  thickness.  These  results  suggest,  however,  that  the  basic  physics  of  the  situation  is 
contained  in  the  equations  and  they  are  sufficiently  accurate  for  many  purposes.  The 
results  are  typical  in  so  far  as  the  wire  resistance  is  comparable  to  the  inductive  reactance 
at  intermediate  frequencies  and  the  wall  thickness  and  conductivity  are  realistic.  These 
are  all  the  comparisons  made  to  date. 

The  code  BOX4  should  predict  the  wire  current  even  in  situations  for  which 
the  analytic  method  is  insufRciently  accurate  to  be  useful.  Because  it  is  comparatively 
inexpensive  to  run  (cf.  section  A.3)  it  may  be  helpful  in  many  problems. 

In  the  remainder  of  section  A.l  the  analytic  method  equations  are  derived  and 
the  basic  characteristics  of  the  wire  current  are  explained  in  both  the  frequency  and  time 
domain.  Sections  A.2-A.5  of  the  report  are  concerned  with  the  description  of  the  code 
BOX4.  The  last  section  discusses  progress  made  in  treating  very  high  frequency  incident 
radiation. 
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A. 2  Diffusion  and  Coupling  Fundamentals 


The  low  frequency  shielding  of  a  metallic  enclosure  can  be  understood  by  con¬ 
sidering  the  case  of  a  spherical  metallic  shell  exposed  to  a  spatially  homogeneous  uniform 
external  magnetic  field  of  intensity  at  frequency  u/  (figure  A-4).  At  frequencies  much 
lower  than  the  natural  modal  frequencies  of  the  sphere  the  displacement  current  can  be 
neglected  (section  A.4)  and  the  internal  field  satisfies 

V  X  H*  =  0  .  (A  -  3) 

Within  the  shell  having  conductivity  a  and  thickness  A  a  current  flows  parallel  to  the 
spherical  surface 

J  =<7E 

where  E  is  the  electric  Integrating  across  th  ...v-r  ness  A  of  the  thin  sheet  yields  the 
surface  current  density  K  (iinh:>eres/meter) 

K  =  oAE  (A  -  4) 

This  current  creates  its  own  magnetic  field  which  subtracts  from  that  of  the  incident 
wave  reducing  the  field  inside  the  shield.  The  change  in  the  tangential  component  of  the 
magnetic  intensity  across  the  shield  is,  of  course 

-  H*  =  K  (A  -  5) 

while  the  component  normal  to  the  shield  is  continuous  through  it.  Finally,  the  electric 
field  in  the  shield  is  associated  with  the  rate  of  change  of  magnetic  flux  in  the  shield  by 
Faraday’s  law 


V  X 


(A -6) 
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It  is  easily  shown  that  the  solution  to  equation  A-3  through  equation  A-6  is  the  internal 
magnetic  held 


H* 


1 

1  _  iwtfngAa 

3 


This  internal  field  is  parallel  to  the  external  field  and  also  spa  tially  homogeneous  (uniform). 
Clearly,  at  frequencies  lower  than 


3 

^iotrAa 


(A -7) 


the  internal  field  almost  equals  the  external  field  because  the  slowly  changing  magnetic 
field  induces  a  small  electric  field  (equation  A-6)  which  drives  little  current  K  in  the  shield. 
But  at  frequencies  higher  than  cji  a  large  current  K  causes  a  small  internal  field.  Above 
the  cut  off  frequency  uii  the  internal  field  decreases  inversely  with  u  and  oscillates  ninety 
degrees  out  of  phase  with  the  external  driver.  This  behavior  is  called  inductive  shielding 
to  distinguish  it  from  the  shielding  which  occurs  at  sufficiently  high  frequencies 


that  the  skin  depth  is  less  than  the  wall  thickness.  In  inductive  shielding  the  electric  field 
through  the  wall  is  nearly  uniform  as  is  the  wall  current  it  drives.  All  but  the  last  section 
of  this  report  is  concerned  with  inductive  shielding  and  frequencies  well  below  both  the 
modal  frequencies  of  the  enclosure  or  wj. 

Several  useful  characteristics  of  the  spherical  case  apply  to  all  enclosures  regard¬ 
less  of  their  shape.  Firstly,  above  a  cutoff  frequency  (approximately  Wi)  the  internal  field 
always  decreases  inversely  with  the  frequency.  Although  the  field  is  not  the  same  at  each 
point  it  decreases  by  the  same  factor  at  each  point  as  the  frequency  increases  (above  the 
cutoff  frequency).  In  other  words,  the  internal  field  is  of  the  form 


H*(f)  =  if(r) 

(jj 

where  f(r)  is  the  spatial  dependence  of  the  field  and  the  frequency  enters  only  as  an  overall 
coefficient.  This  .s  proven  rigorously  for  any  shape  enclosure  in  section  A. 3.  Another 
important  characteristic  is  that  the  internal  fields  are  approximately  spatially  uniform  so 
f(r)  is  approximately  constant.  In  the  cases  of  a  spheie,  cylinder  and  a  parallel  plate 
enclosures  the  internal  field  is  exactly  the  same  at  all  points'*.  This  is  obvious  at  very  low 
frequencies  where  the  walls  do  not  shield  effectively  and  the  internal  and  external  fields  are 
the  same.  As  the  frequency  increases  the  field  in  these  simple  geometries  remains  uniform. 
In  the  general  case,  of  course,  the  fields  will  vary  somewhat  within  the  enclosure  but  as 
long  as  the  shape  is  not  radically  different  from  a  sphere,  cylinder  or  two  parallel  plates 
one  would  expect  the  internal  field  to  be  approximately  uniform.  At  very  high  frequencies, 
of  course,  the  internal  fields  are  very  non-uniform  since  the  cavity  modes  are  excited. 

In  summary,  the  analytic  theory  in  section  A.l  assumes  that  the  field  in  the 
box  is  approximately  spatially  homogeneous  as  it  is  for  a  sphere,  cylinder  and  two  parallel 
plates.  It  might  be  pointed  out  that  the  internal  field  in  these  three  cases  differ  only  by 
about  a  factor  of  two.  The  theory  assumes  that  the  shielding  of  a  box  must  be  quite  similar 
and  the  internal  field  is  approximated  by  the  field  in  a  sphere  having  the  same  volume  as 
the  box  (equation  A-2). 

Equation  A-1  which  estimates  the  wire  current  is  derived  by  integrating  Faraday’s 
law  around  the  dotted  line  closed  path  in  figure  A-2 

/  E  •  d£  =  vjjho  f  H  •  ndS 


The  line  integral  can  be  divided  into  a  part  through  the  wire  and  a  part  through  the  walls 


IR 


/  E-d£-Vd 
Jbca 

where  I  and  R  are  the  wire  current  and  resistance,  respectively.  The  paths  A-B  and  3- 
C-A  are  indicated  in  the  figure.  The  total  internal  magnetic  field  H  is  the  sum  of  the  field 
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produced  by  the  current  in  the  wire  H*^  and  the  field  which  diffuses  through  the  shield  H*. 
Consequently,  H  =  H'*'  +  H*  and 


IR  +  Vd  =  iw/io  jf  IV  •  ndS  +  iuno  J  H*  •  ndS  (A  -  9) 

This  equation  is  proven  rigorously  in  section  A. 2. 

It  was  pointed  out  that  the  magnetic  field  within  the  box  is  roughly  uniform  so 

that 

j  H*  •  ndS  =  H‘Lx 

where  x  and  L  are  defined  in  figure  A-2.  Consider  a  contour  of  integration  in  the  walls 
which  circles  the  box  like  the  wal'  current  (figure  A-2).  The  magnetic  flux  through  the 
contour  is  HwL  so  the  average  field  in  the  walls  must  be 


Ew  = 


/ipwL 
■2(w  +  L) 


iwH- 


Hence, 


Vd  =  Ew(L  +  2x) 


and 


IR  4-  =  iujfio  I H*  •  ndS  4-  iw/zoH'Lx  (A  -  *  0) 

The  flux  produced  by  the  wire  linking  the  path  of  integration  is  proportional  to  the  induc¬ 
tance  £.  of  the  wire 


-Mo/  H*  -ndS 


(A  -  11) 
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Consequently 


IR  + 


■  _  wL(L4-2x)' 

2(w  +  L)  . 


(A  -  12) 


This  is  equation  A-2.  The  authors^  original  derivation  of  this  equation  assumed  that  L 
is  the  inductance  of  a  wire  in  free  space.  Recent  work  has  shown  that  the  free  space 
inductance  is  about  a  factor  of  three  too  large  if  the  wire  is  near  the  center  of  the  box. 
In  section  A.4  the  inductance  is  calculated  accurately.  Figure  A-5  which  graphs  the  in¬ 
ductance  for  the  table  A-2  case  shows  that  it  depends  slightly  on  the  wire  location  and  is 
approximately  the  free  space  inductance.  In  the  section  A.l  comparisons  the  more  accurate 
values  of  inductance  were  employed. 


Clearly,  several  approximations  have  been  made  to  deduce  the  final  results.  The 
most  controversial  are  probably  the  assumption  that  the  internal  fields  are  uniform  and 
approximated  by  a  sphere.  The  authors  believe  these  are  indeed  general  features  of  low 
frequency  diffusion.  The  code  comparisons  support  this  intuitive  judgement  convincingly 
but  much  more  extensive  comparisons  are  needed. 


A. 3  General  Characteristics  of  the  Wire  Current 

Equation  A-9  has  a  .simple  interpretation.  The  right  hand  member  of  the  equa¬ 
tion  is  the  EMF  linking  the  loop  in  figure  A-3.  It  is  the  driver  of  the  wire  current.  The 
wire  inductance  L  is  a  source  of  back  EMF  which  reduces  the  wire  current. The  equation 
indicates  that  the  total  EMF  drops  through  the  wire  resistance  R  and  the  wall  resistance. 
The  path  of  the  loop  in  the  walls  is,  of  course,  arbitrary  but  once  a  path  has  been  chosen 
both  the  driver  and  Vj  must  be  computed  with  that  path. 

The  v/ire  current  dependence  on  frequency  is  typified  by  figure  A-3.  At  low 
frequencies  the  current  is  small  because  the  electric  field  induced  by  the  slowly  changing 
magnetic  field  produces  a  small  EMF  (right  hand  member  of  equation  A-9).  As  the  fre¬ 
quency  increases  the  EMF  and  wire  current  increase.  At  the  frequency  u,  however,  the 
wall  current  begins  to  effectively  reduce  the  magnetic  field  in  the  enclosure  so  the  EMF 
increases  less  rapidly  with  increase  in  frequency.  For  frequencies  larger  than  a  ^  ^ 
the  critical  frequency 
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Inductance  (Henrys  ♦  10^) 


Figure  A-5.  Inductance  as  a  function  of  wire  position  for  table  A  -2  case. 
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u;  »  Wi 


equation  A-2  shows  the  product  of  frequency  and  internal  field  remains  a  constant.  Hence, 
the  EMF  and  wire  current  are  constant  for  frequencies  above  this  value.  Consequently,  the 
maximum  wire  current  occurs  at  frequencies  slightly  above  c^i.  At  much  higher  frequencies, 
of  course,  skin  depth  and  other  effects  alter  this  behavior  and  reduce  the  wire  current  far 
below  its  maximum. 

Equation  A- 12  indicates  that  the  larger  the  wire  resistance  and  inductance  the 
smaller  the  wire  current.  At  low  frequencies 


R 

the  resistance  will  dominate  the  inductive  reactance  and  the  wire  current  will  be  ninety 
degrees  out  of  phase  with  the  internal  magnetic  field.  At  sufficiently  high  frequencies 

R 

w  »  — 

C 

the  resistance  is  negligible  and  the  current  is  in  phase  with  the  internal  field.  As  explained 
above  if  w  <  '<Ji  then  the  internal  magnetic  field  is  in  phase  with  the  external  driving  field 
and  ninety  degrees  out  of  phase  if  a/  2>  C4;i.  The  inductance  can  be  approximated  by  that 
of  a  wire  in  free  space  although  section  A.4  derives  a  more  accurate  value. 

The  wire  current  can  also  be  derived  in  the  time  domain.  Substituting  equa¬ 
tion  A-1  in  equation  A-2  yields 


I(a;) 


iwGH*^  (w) 

(R  —  iuj£){l  —  iu/uji) 


where 


G  —  flQ 


xL 


wL(L  +  2x) 
2(w  +  L) 
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This  is  the  wire  current  as  a  function  of  the  external,  driving  magnetic  field 
in  the  frequency  domain.  By  taking  the  inverse  fourier  transform  the  wire  current  I(t)  as 
a  function  of  time  is  found  to  be 


This  is  the  wire  current  as  a  function  of  time  driven  by  an  arbitrary  external  pulse  whose 
time  dependence  is  (t). 

Clearly,  the  wire  current  depends  upon  Wi  and  the  characteristic  frequency  R/il. 
The  first  is  the  cutoff  frequency  below  which  the  enclosure  shields  the  external  field  very 
poorly.  The  second  is  the  time  constant  of  the  R  -  circuit  formed  by  the  wire  and  the 
box.  Fluctuations  faster  than  cji  and  R/£  are  smoothed  out  by  this  integral.  In  general, 
when  the  driver  is  first  turned  on  the  current  slowly  builds  in  the  wire  and  after  it  is  turned 
off  the  wire  current  slowly  decays  to  zero. 

The  authors  believe  that  the  frequency  domain  solution  in  equations  A-1,  and  A- 
2  and  the  above  time  domain  solution  contain  the  basic  physics  of  the  problem  and  permit 
reasonably  accurate  predictions  in  most  situations.  A  straightfoward  extension  of  the 
theory  will  enable  predictions  to  be  made  for  wires  of  any  shape. 


A.4  RESULTS  AND  BASIC  PHYSICS 


Figure  A-2  illustrates  a  metallic  enclosure  driven  by  an  external  magnetic  field 
directed  along  one  of  its  edges  and  oscillating  at  the  frequency  u.  Because  the  wavelength 
is  much  longer  than  the  box  dimensions  the  field  is  independent  of  position.  If  the  current 
on  the  internal  wire  can  be  predicted  for  this  case  then  it  can  be  predicted  for  any  long 
wavelength  external  driving  field.  It  is  only  necessary  to  fourier  analyze  thn  more  complex 
driver  into  its  frequencies  and  resolve  the  field  at  each  frequency  into  its  components  along 
the  three  orthogonal  edges.  The  total  wire  current  is  the  sum  of  the  currents  produced  by 
each  part  of  the  driver. 

The  electromagnetic  field  which  is  produced  by  the  current  within  the  wire  obeys 
Faraday’s  and  Ampere’s  laws 


V  X  E’*  =  iw/ioH* 
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V  X  H'^  =  J*  -  iwcoE 


within  the  box  where  J*'  is  the  current  density  in  the  wire.  A  single  frequency  is  considered 
which, as  explained  above,  is  not  a  loss  of  generality.  In  this  report  all  quantities  are 
assumed  to  vary  as  exp(-i  <jj  t).  Within  the  walls  having  conductivity  a 


V  X  E*'  = 


V  X  H*  =  oE*'  -  iwcoE* 


and  outside  the  enclosure  the  fields  satisfy  both 


V  X  H'^  =  -iweoE* 


(A  - 13) 


and  Faraday’s  law. 

The  box  is  radiated  with  a  driving  field  which,  presumably,  is  produced  by  an 
external  current  density  J*.  The  current  density  might,  for  instance,  be  the  currents  within 
a  HPM  microwave  weapon  or  the  Compton  current  in  the  atmosphere  which  generates  an 
EMP  pulse.  J*  is  assumed  to  be  known.  Outside  the  box  the  field  produced  by  this  driver 
satisfies  Faraday’s  law  and 


V  X  H*  =  J*  -  iwcoE*  . 
In  the  walls  the  governing  equation  is 

V  X  H*  =  oE*  -  iwcoE* 


and  within  the  box 


V  X  H*  =  -iweoE* 


(A  - 14) 
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The  presence  of  ^he  wire  is  correctly  omitted  in  equation  A- 14 

The  total  electric  Et  and  magnetic  Hr  fields  are  the  sum  of  the  fields  produced 
by  the  external  driver  J*,  the  current  in  the  wire  J*'  and  the  current  in  the  walls  aEr 


Et  =  E*  +  E‘ 


Ht  =  H’  +  H* 


However,  the  wire  cui-.dnt  J*'  is  caused  by  the  total  electric  field  according  to 


J*  =  <7''(E*'  +  E‘)  (A  -  15) 

in  terms  of  the  wire  conductivity,  <7*.  Solving  equations  A- 13  through  A- 15  furnishes  the 
current  J’^in  the  wire.  The  wire  current  derived  in  tnis  way  is  exact  but  it  is  difficult  to 
solve  the  equatiorxs  if  the  wire  is  so  thick  tin  t  the  current  around  it  is  non-uniform.  The 
authors  will  attempt  this  less  important  case  iu  the  future. 

This  approaych  to  the  problem  which  was  pointed  out  by  Dr  James  Gilbert  at 
Mission  Research  Corporation  is  similar  to  the  technique  used  to  solve  many  scattering 
problems.  Solving  equations  A-13  yields  the  electromagnetic  fields  produced  oy  a  specified 
wire  current  J*  when  the  external  driver  J*  is  not  present.  Solving  equations  A-14  yields 
the  electiomagnetic  fields  due  to  the  known  external  diiver  J'  when  the  wire  is  not  present. 
These  two  aiinplei  problems  are  solved  in  sections  A. 3  and  A. 4.  Equation  A-15  then  yields 
the  wire  cuirrent  trivially.  Of  course,  in  equation  A-15  E*  is  a  function  of  the  unknown 
J”.  So  2-3  is  actually  an  implicit  equation  which  must  be  solved  to  find  J*'.  The  authors 
employed  this  method  to  derive  the  wire  current  results  in  the  previous  section. 

The  electric  field  can  be  written  in  terms  of  the  scalar  and  vector  potentials^,  (p 

and  A 


E’'  =  iwA*  -  V(^* 

E*  =  iwA*  -  (A  -  16) 
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Integrating  equation  A- 15  along  tn#*  wire  (figure  A-2)  yields 


IR 


=  E*  +  E*dx  =  <^1  -  <f>^  +  ~  (i>B  +  r  A^'dx  -}-  iu  A*dx  (A  -  17) 

^  A  » A  'A 


where  R  is  the  wire  resistance. 

Consider  the  integration  path  through  the  walls  shown  in  figure  A-2.  The  line 
integral  of  the  electric  field  along  this  path  is  called 


Vd  =  y  Ejd^  =  +  iu  j  A*d£  +  iw  ^  A* 

adding  equation  A-17  and  equation  A-18  yields 


d£ 


(A  -  18) 


IR  +  Vd  =  iw  /  A*d<  +  iu  /  A’d<  (A  -  19) 

Where  the  closed  path  is  formed  by  the  wire  and  the  path  through  the  walls.  But  the  line 
integrals  can  be  related  to  the  surface  integral  of  the  magnetic  field  through  the  closed 
path 


y  A  •  d£  =  /{ V  X  A)  •  ndS  =  ^  B  •  ndS 
by  Stokes  theorem  yielding 

IR  +  Vd  =  iw  y  B*  •  ndS  +  iij  Jb'-  ndS  (A  -  20) 

This  is  a  rigorous  proof  of  equation  A-9,  the  starting  point  for  the  approximate  wire  treat¬ 
ment  in  section  A.l.  It  is  important  to  note  that  in  adding  equation  A-17  and  equation  A-18 
the  scalar  potentials  drop  out.  Consequently,  equation  A-20  is  correct  even  in  the  general 
case  for  which  the  scalar  potential  is  not  negligible.  Calculating  the  wire  current  with  this 
equation  requires  that  Vd  be  estimated,  however,  a  complication  which  is  avoided  in  the 
following  manner. 

The  scalar  potentials  are  small  and  can  be  neglected  in  equation  A-17.  This 
is  shown  in  sections  A. 3  and  A.4  and  explained  intuitively  in  the  following  way.  In  the 
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Lorentz  gauge  the  scalar  and  vector  potentials  are  driven  by  the  charge.^  and  current 
density  J,  respectively^ 


vV 


1  av 

c*  at* 


1 


— p 

Co 


V*A- 


1  a*A 
c*  at* 


—  —f^oJ 


In  general,  charge  accumulates  on  the  surface  of  the  wire  and  creates  an  scalar  potential 
whose  gradient  is  the  irrotational  part  of  the  electric  field.  The  charge  accumulates  in  such 
a  way  that  the  boundary  condition  (equation  A-15)  is  satisfied  on  the  wire.  For  instance,  if 
the  wire  is  coiled  on  a  cylinder  to  form  an  inductor  then  a  small  charge  density  is  deposited 
by  the  almost  divergenceless  wire  current  on  the  surface.  For  a  zero  resistance  wire  the 
field  on  its  surface  must  be  zero  so  the  electric  field  ituA  must  be  equal  and  opposite  to  that 
produced  by  the  scalar  potential  V<^.  In  a  small  number  of  very  symmetrical  situations, 
however,  the  vector  potential  alone  is  sufficient  to  satisfy  the  boundary  condition. 

Consider,  for  example,  a  zero  resistance  wire  loop  in  a  homogeneous  oscillating 
magnetic  field  (figure  A-6)  which  induces  an  EMF  around  the  wire  and  causes  a  current 
to  flow.  The  oscillating  current  creates  an  electric  field  iwA  at  the  wire  which  just  cancels 
the  induced  field.  It  is  easily  shown  that  charge  does  not  accumulate  along  the  thin  wire 
and  the  scalar  potential  is  zero.  But  if  the  circular  loop  is  distorted  slightly  then  the 
vector  potential  can  not  by  itself  zero  the  electric  field  and  a  scalar  potential  must  be 
produced.  Similarly,  if  the  oscillating  driving  field  is  not  homogeneous  the  scalar  potential 
is  significant. 

The  vector  potential  A*  is  the  sum  of  the  potentials  produced  by  the  current  in 
the  walls  A'  and  the  external  source  A*^.  The  low  frequency  external  source  produces  a 
homogeneous  magnetic  field  directed  along  the  z  axis  (figure  A-2)  so  that 


Hd  =  —  V  X  A° 
tM) 

Since  the  vector  potential  is  symmetrical  about  an  socis  through  the  center  of  the  box 
parallel  to  the  z  axis  the  correct  solution  to  this  equation  is 
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H 


Figure  A>6.  Wire  in  an  oscillating  magnetic  Held. 

Aj  =  ^-(y->o 

=  (A -21) 

Hence,  equation  A- 17  becomes 

lR  =  \uf  j  A*dx  +  \ujJ  A*dx  +  iu  j  A°dx  (A  —  22) 

This  final  equation,  rather  than  equation  A-17,  is  used  to  find  the  wire  current  with  the 
code. 


86 


In  summary,  the  method  which  the  authors  developed  to  calculate  accurately 
the  wire  current  er.ploys  the  computer  code  discussed  in  section  A. 3  to  calculate 


A'^dx 


and  the  analytic  calculation  in  section  A.4  to  yield 


Then  equation  A-21  yields 


/V'-dx 

Since  the  magnetic  field  is  proportional  to  the  wire  current  I 

A'"dx  =  £1  (A  -  23) 

Finally,  the  wire  current  is  obtained  by  solving  the  algebraic  equation  A-22  which  becomes 

IR  +  \<jjC1  =  i(jJ  J  A'dx  +  iw  A’^dx  (A  —  24) 

It  is  worthwhile  pointing  out  that  C  is  the  inductance  of  the  wire  in  the  enclosure 
in  the  limit  that  the  walls  shield  the  outside  of  the  box  from  the  wire  very  effectively.  This 
equation  differs  from  equation  A-20  which  is  the  starting  point  of  the  analytic  method  in 
two  ways.  Firstly,  it  is  not  necessary  to  calculate  Vd,  the  line  integral  of  the  electric  field 
in  the  walls.  Secondly,  the  scalar  potential  has  been  neglected  in  equation  A-23  but  it  is 
included  in  the  surface  integrals  in  equation  A-20. 


A.5  EXTERNAL  DRIVER  CALCULATION 


In  this  section  the  electric  field  produced  at  the  wire  by  the  external  magnetic 
field  which  diffuses  through  the  shield  is  derived.  Outside  the  wall  the  magnetic  intensity 
H  can  be  expressed^  by  the  gradient  of  a  magnetic  scalar  potential  Cl 
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H  =  vn 


where 


V*n=:0  . 


It  is  well  known^  that  within  the  enclosure  the  scalar  potential  can  be  expressed 
as  a  function  of  its  value  and  derivative  on  the  inside  surface  of  the  wall 


(A  -  25) 


where 


G  = 


1  1 
4n  |r  —  r' 


Similarly,  outside  the  rectangle 


n~>(r)  =  -  /  (G^  -  ||n~*)ds'  (A  -  26) 

where  the  normal  derivative  is  always  outside  the  box  and  is  the  external  field  imposed 
on  the  metallic  enclosure  by  the  source  such  as  a  low  frequency  electromagnetic  wave. 

The  current  in  the  thin  walls  is  expressed  as  a  surface  current  density  K  (am¬ 
peres/meter)  parallel  to  the  conducting  surface.  Assuming  it  is  divergenceless  permits  it 
to  be  expressed  in  terms  of  a  scalar  stream  function  ^ 


K  =  n  X  V.V'  (A  -  27) 

where  n  is  the  surface  normal  and  the  surface  divergence  is  V,. 

Latham  and  Lee^  have  shown  the  magnetic  scalar  potentials  on  the  inside  and 
outside  surfaces  of  the  walls  obey 
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f  /  GV;VdS' 

J  on!  iuiuoZ  J 


iWMog 


-n°“*  =  n“  +  /  n°“*-~ds'  -  — —  f  Gv;Vds'  (a  -  28) 

2  J  dn>  lUfiQg  y  *  ^  ^ 

where  g  =  a  A,  the  product  of  wall  conductivity  and  thickness. 

Making  use  of  the  facts  that  the  normal  derivative  of  H  is  continuous  across  the 
wall  whereas  the  jump  in  the  tangential  component  of  H  across  the  wall  is  equal  to  the 
surface  current  density  K  yields  an  equation  for  the  stream  function* 


^G  dG 
dn"  dn' 


dS'dS" 


-  -  '^(0^;GdS'  (A  -  29) 

lufuog  J 

This  integral  equation  is  ideally  suited  to  numerical  solution  since  only  integrals  over  the 
unknown  tp  rather  than  derivatives  appear.  The  function  tp  is  calculated  at  a  several 
hundred  locations  on  the  surface  of  the  box.  For  the  i***  location  the  integral  equation  is 
approximated  by 


4^  4^  dn' 


+  EE-/'k 

;  k 


5ni,  dnj 


1 

iw^tog 


where  £*  is  the  area  of  each  surface  cell.  This  yields  a  system  of  several  hundred  simulta¬ 
neous  linear  algebraic  equations  which  are  then  solved  numerically  with  the  code  BOX4 
io  yield  the  stream  function.  Once,  xp  is  known  the  magnetic  scalar  potential  H  and  the 
magnetic  field  can  be  calculated  within  the  box  with  little  additional  difficulty.  Further- 
Uiore,  the  magnetic  vector  potential  A  at  the  wire  produced  by  the  current  in  the  walb 
can  be  found  by  integrating  over  the  stream  function 
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(A  -  30) 


—  Hq  j  dS'G{n  X  V,0) 

This  is  the  quantity  needed  in  section  A.4  (equation  A-24)  to  calculate  the  wire  current. 

The  numerical  solution  is  fast  and  reasonably  accurate.  In  the  high  frequency 
limit  the  walls  shield  the  external  magnetic  field  very  well.  Consequently,  the  magnetic 
vector  potential  A^  (section  A.4)  produced  by  the  known  driver  is  equal  and  opposite  to 
that  produced  by  the  wall  current  (equation  A-30).  In  the  figure  A-3  example  they  differed 
by  only  eight  percent  when  0  was  approximated  by  96  ceils.  This  accuracy  was  achieved 
for  a  total  computer  cost  of  about  fifty  cents.  Dividing  the  box  into  216  cells  achieved 
slightly  better  accuracy  for  a  cost  of  about  ten  dollars.  The  computations  were  done 
on  an  ELXSI  time-sharing  computer  which  is  similar  to  the  popular  Digital  Equipment 
Corporation  VAX780.  Most  of  the  computer  time  is  spent  evaluating  the  coefficients  in  the 
simultaneous  equations  many  of  which  are  multiple  integrals  rather  than  actually  solving 
the  equations  which  is  done  by  the  Gaussian  elimination  method.  A  significant  savings  in 
computer  cost  could  be  realized  by  making  better  use  of  the  problem’s  symmetry  when 
computing  the  multiple  integrals. 

An  important  simplification  is  understood  by  examining  the  limit  of  equation  A- 
28  as  the  frequency  becomes  much  larger  than  wi.  Then  the  last  terms  in  these  equations 
are  small  and  to  zero  order  the  internal  field  D'"  is  zero  since  the  walls  shield  well  at  high 
frequencies.  Equation  A-28b  yields  a  zero  order  approximation  to  the  external  potential 
Qout  can  be  substituted  back  into  equation  A-28a  to  yield  a  better  estimate  for  the 

internal  field 


-n‘'‘(f)  +  f  n‘"(r)^ds'  = .  ---  f  Gv*n»“‘ds' 

2  ^  ^  J  ^  '  dn  iuJUoaA  J 

The  scalar  potential  IT"  is  a  function  of  frequency  and  position  r.  But  the  frequency 
appears  only  in  the  coefficient  of  the  right  hand  member.  Consequently,  as  the  frequency 
is  increased  the  field  H'"  at  each  location  decreases  inversely  with  the  frequency.  Hence,  * 
the  vector  potential  A*  decreases  at  the  wire  inversely  with  frequency.  As  was  pointed  out 
in  section  A.l  this  has  important  consequences. 

The  assumption  that  the  wall  current  is  divergenceless  implies  that  there  is  no 
charge  accumulation  on  the  walls.  Therefore,  the  scalar  potential  <l>  must  (under  this 
assumption)  be  zero.  As  explained  in  section  A.4,  usually,  this  is  not  true.  In  this  approach 
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it  is  impossible  to  calculate  the  contribution  to  the  electric  field  produced  by  the  scalar 
potential.  Fortunately,  BOX4  shows  that  the  vector  potential  A‘  produced  by  the  wall 
current  is  approximately  constant  along  the  wire.  Section  A.6  indicates  that  the  potential 
produced  by  the  wire  A*'  is  also  constant.  But  the  total  electric  held  on  the  wire  must  be 
constant  since  it  equals  IR.  Since  A*^  and  A*  are  constant  they  are  sufficient  to  satisfy  the 
boundary  condition  on  the  wire  and  <}>  must  be  negligibly  small. 

Toward  the  edges  of  the  wire  A*^  is  not  constant.  The  authors  believe  this  is  due 
to  numerical  inaccuracy.  Further  attention  to  this  point  is  necessary. 


A.6  WIRE  CURRENT 


As  the  oscillating  current  flows  through  the  wire  and  around  the  walls  of  the 
enclosure  it  produces  a  substantial  electric  held.  If  the  wire  resistance  is  zero,  the  mag¬ 
nitude  of  the  field  (along  the  wire)  which  penetrates  the  enclosure  is  equal  to  that  of  the 
field  produced  by  the  wire  current.  In  this  section  the  fields  produced  in  the  box  by  the 
wire  current  are  calculated  so  that  the  integral  in  equation  A-23  and  the  inductance  can 
be  evaluated. 

Figure  A-1  illustrates  a  rectangular  enclosure  driven  by  a  current  I  which  runs 
through  its  central  wire  and  returns  via  the  metallic  walls.  The  current  oscillates  at  the 
frequency  w  at  a  constant  amplitude.  It  is  not  necessary  to  specify  the  source  of  energy 
which  drives  the  current  against  the  ohmic  losses  in  the  walls  and  the  wire.  In  the  limit  of 
very  low  frequency,  small  conductivity  or  small  wall  thickness  the  current  will  distribute 
itself  so  that  the  electric  field  at  the  walls  is  sufficient  to  drive  the  current.  In  general, 
appreciable  magnetic  fields  will  leak  outside  the  enclosure.  But  as  the  conductivity  or  wall 
thickness  increase  the  electric  field  at  the  walls  decrease.  Furthermore,  as  the  frequency 
increases  the  electric  field  induced  within  the  enclosure  by  the  changing  current  becomes 
much  larger  than  that  necessary  to  drive  the  wall  current.  In  either  limit  the  boundary 
condition  at  the  walls  is  that  the  electric  field  at  the  walls  becomes  negligible.  Consequently, 
the  wall  conductivity  can  be  considered  infinite.  The  conditions  for  which  this  idealization 
is  valid  are  quantified  at  the  end  of  this  section.  In  the  case  of  the  conductivity  and 
wall  thickness  referred  to  in  the  section  A.l  examples  the  approximation  is  justified  for 
frequencies  above  about  one  kilohertz. 
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The  vector  potential  within  the  infinite  conductivity  box  satisfies 


V*A  +  -jA  =  -/ioJ 


(A  -  31) 


which  has  the  solution 


X  =  j  G(f,r/)JdV 

where  the  Green’s  function  satisfies  the  appropriate  boundary  conditions  on  the  walls.  J  is 
the  wire  current  density.  In  terms  of  the  eigenfunctions  satisfying  the  Helmholtz  equation 

W.  +  =  0  (A  -  32) 


the  Green’s  function^  is 


G  =  Moc’E 

n 


w’  -  w* 


(A  -  33) 


Equation  A-33  indicates  that  if  the  frequency  u*  is  much  less  than  the  natural 
modal  frequencies  of  the  box  u)n  then  the  vector  potential  does  not  differ  appreciably  from 
the  value  if  would  have  if  a;  =  0.  Consequently,  it  is  only  necessary  to  solve  the  quasi-static 
equation 


V*A  =  -noJ 


subject  to  the  boundary  conditions  that  the  normal  component  of  the  magnetic  field  on 
the  wall  is  zero.  The  eigenfunctions  are  then 


0inn 


miry.  .  .nxz. 
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and  the  vector  potential  is 


A* 


8/^1 

bcjT* 


E 


sin(2ft)8in(2p) 

(i)’ + (?)» 


8in{-j— )8m(— )Qm„ 

D  C 


(A  -  34) 


where  the  wire  position  is  y  =  ^  and  z  =  ri  and 

Q»  =  «p{-(f)»[(f)'  +  (f)']}  . 

It  is  important  to  consider  a  wire  of  finite  radius  to  calculate  the  wire  induc> 
tance.  Furthermore,  the  series  converges  very  poorly  if  the  wire  radius  is  zero  since  this 
necessitates  very  high  order  spatial  modes  because  the  vector  potential  changes  direction 
across  the  wire.  In  this  report  the  wire  current  density  J  has  been  assumed  to  fall  off  as  a 
Gaussian  with  characteristic  wire  radius  wo 


From  the  above  series  expression  for  A  the  magnetic  flux  through  the  contour  in 
figure  A-2  can  be  derived 


lev  X  A  •  n)dS 


It  is  the  negative  of  the  electromotive  force  produced  around  the  loop  by  the  wire  current. 
Since  the  walls  are  infinitely  conducting  it  is  also  equal  to 


EMF  = 


iw  ^  A  •  dx 


(A  -  36) 


These  are  the  final  expressions  needed  in  section  A.4  (equation  A-23)  to  calculate  the  wire 
current.  The  inductance  £  is  defined  in  terms  of  this  integral  by  equation  A-23. 

The  circumstances  for  which  the  infinite  conductivity  approximation  is  valid  are 
easily  quantified.  The  resistance  of  a  conducting  sheet  whose  length  and  width  are  equal 
is 
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R, 


1 

(tA 


(A  -  36) 


where  a  and  A  are  the  conductivity  and  thickness.  For  instance,  a  sheet  of  conductivity 
10^  mho/m  and  one  millimeter  thickness  has  a  resistance  of  0.001  ohms.  Clearly,  if  an  ohm 
meter  is  used  to  measure  the  resistance  of  the  walls  between  points  A  and  B  (on  figure  A-2) 
a  resistance  of  this  order  of  magnitude  will  be  observed  provided  the  wall  conductivity  and 
thickness  are  the  same  as  the  sheet  and  the  box’s  length,  width  and  height  are  comparable. 
The  line  integral  of  the  electric  field  in  the  wall  from  point  A  to  B  is  path  dependent  but 
is  typically  of  the  order  IRw  in  terms  of  the  wire  current  I.  The  EMF  around  the  entire 
loop  in  figure  A*^2  is  shown  in  expression  A-3S.  The  infinite  conductivity  approximation  is 
valid  if  the  EMF  greatly  exceeds  IRw  or 

IRw  <  u;  A  •  dx  (A  -  37) 

It  has  been  assumed  in  equation  A-35  that  the  contribution  to  the  electric  field 
due  to  the  scalar  potential  is  negligible.  It  is  easily  seen  from  equation  A-35  that 

V-A  =  0 

Consequently,  from  the  Lorentz  condition  the  scalar  potential  is  also  zero.  Therefore, 
for  frequencies  much  lower  than  the  modal  frequencies  of  the  enclosure  (frequently  below 
hundreds  of  megahertz)  but  sufficiently  high  that  the  walls  provide  reasonable  magnetic 
field  shielding  (frequently  above  one  kilohertz)  the  scalar  potential  causes  little  electric 
field.  Hence,  there  is  no  charge  in  the  system  and  the  flow  of  current  is  divergenceless. 
It  is  important  to  understand  that  at  higher  frequencies  the  scalar  potential  can  not  be 
neglected.  Furthermore,  the  straight  wire  is  a  particularly  simple  case.  If  the  wire  were 
to  lead  to  a  circuit  containing  inductors  or  capacitors  then  these  would  produce  a  scidar 
potential.  If  these  circuit  components  occupied  a  large  fraction  of  the  volume  of  the  box 
the  computer  approach  taken  in  this  report  would  require  modification. 

The  inductance  has  not  been  carefully  computed  at  low  frequencies  whre  equa¬ 
tion  A-37  is  not  valid. 
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A.7  CALCULATION  OF  HIGH  FREQUENCY  RESPONSE 


In  the  previous  sections  we  have  described  in  detail  our  w^ork  on  computation  of 
the  low  frequency  response  of  the  wire,  in  particular  the  solution  of  the  interior  potential 
equation,  including  the  wire.  In  this  section  we  shall  discuss  the  problem  of  high  frequency 
response  calculations.  Since  the  low  frequency  solution  was  the  main  thrust  of  our  early 
efforts  the  progress  here  has  been  more  limited.  To  date  we  have  almost  completed  the 
development  of  the  code. 

The  surface  current  on  a  box  illuminated  by  electromagnetic  radiation  can  be 
described*  in  terms  of  the  incident  magnetic  fields.  A  careful  analysis*,  yields  the  equation: 

K{i)-2lnxKx  G{r,f')d8'  =  2H“‘{f)  (A  -  38) 

where  K  is  the  surface  current,  is  the  external  magnetic  field  and  G  is  the  outgoing 
Green  function 


G(r,f)  =  1/  {47r|f- r|  exp{ik|r  -  ?!)}  (A  -  39) 

corresponding  to  frequency  w  =  k  c. 

The  above  equations  provide  the  current  density  on  the  surface  of  the  box.  Ref¬ 
erence  1  indicates  how  these  are  used  to  drive  the  internal  fields  and  wire  current. 

In  our  first  effort  we  have  focussed  on  the  solution  of  equation  A-38,  i.e.  the 
determination  of  the  surface  current  density  on  the  outside  of  the  box.  Before  we  proceed 
to  discussion  of  the  numerical  solution  of  the  problem  we  note  that  ‘‘high  frequency”  in  our 
case  usually  indicates  frequencies  of  1  MHz  and  higher.  Since  the  integral  in  equation  A-38 
is  oscillatory  in  nature,  the  required  number  ofsteps  for  the  calculations  is  dictated  by  the 

*  standard  assumption  of  10  -  20  points  per  period. 

To  solve  the  problem  numerically  we  have  adopted  the  standard  gridding  of  the 

*  box  faces  as  indicated  in  our  solution  of  the  low  frequency  problem.  For  a  system  of 
equations  to  have  a  meaningful  result  typical  matrix  dimensions  should  not  exceed  1000 
X  1000  or  so.  This  limits  the  number  of  cells  on  each  face  to  42  or  about  6  cells  in  each 
direction.  This  implies  that  our  solution  will  be  quite  accurate  for  a  box  whose  size  is 
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smaller  than  about  third  of  a  wavelength.  For  the  1  MHz  frequency,  this  corresponds  to 
boxes  of  about  100  meters  across. 

For  larger  boxes,  the  solution  will  naturally  be  less  accurate  and  we  shall  use  our 
code  to  provide  an  initial  es  Jmate  of  the  surface  current.  This  solution  will  be  improved 
upon  by  inserting  it  into  a  very  accurate  integrator  of  the  LHS  of  equation  A-38  and  com-  « 
paring  with  the  RHS  driver.  Such  integrator  is  not  constrained  by  memory  requirements 
and  does  not  include  any  matrix  operations. 

We  expect  these  two  codes  will  be  sufficient  to  provide  a  general  trend  prediction 
for  the  solutions  computed  but  we  are  aware  that  the  solution  accuracy  is  going  to  be 
limited.  However,  since  our  goal  is  to  develop  simple  models  of  surface  currents,  we  feel  that 
the  finite  difference  approach  should  be  adequate.  Furthermore  the  accurate  integration  of 
the  LHS  of  equation  A-38,  once  coded,  will  provide  a  test  method  for  evaluation  of  models 
of  the  surface  current  density. 

After  expansion  of  the  vector  product  and  if  we  replace  the  integral  by  a  sum, 
the  integral  equation  becomes  a  matrix  equation  for  the  surface  currents.  For  the  present 
problem  we  have  considered  only  a  rectangular  box.  Other  box  geometries  could  be  also 
computed  but  the  algebra  becomes  somewhat  more  involved.  For  the  six  faces  of  the  box, 
there  are  12  independent  components  of  the  current  together  with  the  boundary  conditions 
at  the  edges  of  the  box  which  arise  from  continuity.  Each  component  is  complex  and  the 
physical  current  is  the  real  part  of  the  product  K  exp(iwt). 

The  resultant  matrix  equation  can  be  cast  in  the  form; 

A1=C 

where  A  is  the  LHS  matrix  (kernel)  and  D  is  the  driving  term.  In  order  to  save  memory, 
we  had  hoped  that  the  kernel  has  some  inherent  symmetries  that  permit  more  efficient 
computation.  However,  as  it  turns  cut,  the  matrix  A  is  band  structured,  but  not  symmetric 
(or  antisymmetric).  Thus  the  computation  of  the  surface  currents  requires  large  memory  ^ 
to  be  available. 

Currently,  the  coding  is  essentially  completed,  with  some  problems  occurring 
with  the  mixing  of  real  and  complex  arithmetics.  We  expect  to  have  some  results  within  * 
the  next  couple  of  weeks.  We  shall  spend  a  significant  amount  of  time  verifying  the  code 
results  for  the  limiting  cases. 


96 


The  next  stage  of  our  calculations,  solving  the  Maxwell  equations  in  the  box 
interior  can  be  accomplished  with  codes  already  available.  This  effort  is  expected  to  be 
relatively  straightforward. 

Our  next  task  is  the  development  of  approximate  formulas  for  the  surface  cur¬ 
rents.  While  it  is  possible  to  rely  solely  on  numerical  experiments  we  are  also  conducting 
an  analysis  of  the  structure  of  the  solution  and  a  survey  of  the  literature  for  available  exact 
solution  in  special  geometries. 
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The  approach  involved  frequency  domain  analysis  based  on  three  distina  frequency  regimes: 
low,  intermediate,  and  high.  In  the  low  frequency  region  the  skin  depth  is  larger  than  the  wall  thick¬ 
ness  of  the  enclosure.  The  intermediate  region  extends  to  the  lowest  internal  mode  of  the  enclosure. 
The  high  frequency  region  extends  upward  from  the  lowest  mode.  Validation  of  the  analytical  results, 
planned  under  Phase  II  of  the  SBIR,  were  not  performed. 
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